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Abstract 

We consider an Allen-Cahn type equation of the form ut = Au + 
£~^/^(x, t, u), where £ is a small parameter and f'^{x,t,u) = f(u) — 
eg^{x^t,u) a bistable nonlinearity associated with a double- well po- 
tential whose well-depths can be slightly unbalanced. Given a rather 
general initial data ug that is independent of e, we perform a rigor- 
ous analysis of both the generation and the motion of interface. More 
precisely we show that the solution develops a steep transition layer 
within the time scale of order e^|lne|, and that the layer obeys the 
law of motion that coincides with the formal asymptotic limit within 
an error margin of order e. This is an optimal estimate that has not 
been known before for solutions with general initial data, even in the 
case where g'^ = 0. 

Next wG consider systems of reaction-diffusion equations of the form 

Ut = Am + e"^ /"^(w, v) 
Vt = DAv + h{u, u), 

which include the FitzHugh-Nagumo system as a special case. Given 
a rather general initial data (uo,wo), we show that the component u 
develops a steep transition layer and that all the above-mentioned re- 
sults remain true for the M-component of these systems. 

Key Words: nonlinear PDE, reaction-dilfusion system, singular pertur- 
bation, Allen-Cahn, FitzHugh-Nagumo, interface motion 0. 
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1 Introduction 



1.1 Perturbed Allen-Cahn equation 

In some classes of nonlinear diffusion equations, solutions often develop in- 
ternal transition layers — or "interfaces" — that separate the spatial domain 
into different phase regions. This happens, in particular, when the diffusion 
coefficient is very small or the reaction term is very large. The motion of 
such interfaces is often driven by their curvature. A typical example is the 
Allen-Cahn equation ut = Au + e~'^ f{u), where e > is a small parameter 
and f{u) is a bistable nonlinearity, whose meaning will be explained below. 
A usual strategy for studying such phenomena is to first derive the "sharp 
interface limit" as e — > by a formal analysis, then to check if this limit 
gives good approximation of the behavior of actual layers. 

In this paper we study a perturbed Allen-Cahn type equation of the form 



^ Ut = Azi + ^(/(u) — eg^{x,t,u)) in 17 x (0, +cx)) 

Fill 

— = on X (0, +CX)) 

ov 

^u{x^Q) = uq{x) in $7, 



and study the behavior of layers near the sharp interface limit as e — > 0. Here 
is a smooth bounded domain in (A^ > 2) and u is the Euclidian unit 
normal vector exterior to d^. The nonlinearity is given by f[u) := — W'{u), 
where W{u) is a double-well potential with equal well-depth, taking its 
global minimum value at n = a±. More precisely we assume that / is 
and has exactly three zeros a_ < a < a-)_ such that 

f\a±) < 0, /'(a) > (bistable nonlinearity), (1-1) 

and that 

f{u)du = 0. (1.2) 



The condition (jl.ip implies that the potential W{u) attains its local minima 
at u = Q;_,a+, and (jl.2p implies that W{a-) = W{a+). In other words, 
the two stable zeros of /, namely a_ and a+, have "balanced" stability. A 
typical example is the cubic nonlinearity /(n) = n(l — n^). 

The term eg^ represents a small perturbation, where g^{x, t, u) is a func- 
tion defined on x [0, +oo) x M. This has the role of breaking the balance 
of the two stable zeros slightly. In the special case where g^ = 0, problem 
(P^) reduces to the usual Allen-Cahn equation. As we will explain later, 
our main results are new even for this special case. 

We assume that g^ is in x and in t, u, and that, for any T > 
there exist ?9 G (0, 1) and C > such that, for ah {x, t,u) gIIx [0, T] x M, 

\A^g%x,t,u)\ <Ce~^ and \gf{x,t,u)\<Ce-^, (1.3) 
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\gl{x,t,u)\<C, (1.4) 
||/(-,-,n)|| 1+^ _ < C. (1.5) 

Moreover, we assume that there exists a function g{x^ t, u) and a constant, 
which we denote again by C, such that 

|/(a;, t, u) - g(x, t,u)\ < Ce, (1.6) 



for ah smah e > 0. Note that the estimate (jl.Sp and the pointwise conver- 
gence ^ g (as e ^ 0) imply that g satisfies the same estimate as (|1.5|) . 
For technical reasons we also assume that 



^ = ondnx [0, Tl X M, (1.7) 
du 

which, in turn, implies the same boundary condition for g. Apart from these 
bounds and regularity requirements, we do not make any specific assump- 
tions on the perturbation term g^. 

Remark 1.1. Since we will consider only bounded solutions in this paper, it is 
sufficient to assume (ll.3p - (ll.5p to hold in some bounded interval —M < u < 
M. Note that if g^ does not depend on e, then the assumptions (|1.3p - (|1.5p 
are automatically satisfied on any interval —M<u<M. □ 

Remark 1.2. The reason why we do not assume more smoothness on g is that 
we will later apply our results to systems of equations including FitzHugh- 
Nagumo system, in which g^ loses C^'-'^-smoothness as e — > 0. □ 

As for the initial data uo{x), we assume uq S C'^{Q). Throughout the 
present paper the constant Cq will stand for the following quantity: 

Co := \\uo\\co(n) + l|Vno||cO(n) + ll^^^ollco(n)- (1-8) 

Furthermore we define the "initial interface" Fq by 

Fo := {x e $7, Uo{x) = a}, (1.9) 

and suppose that Fq is a C^^^ hypersurface without boundary such that, n 
being the outward unit normal vector to Fq, 

FoCCO and Vno(a;) • n(x) / ifxGFo, (1.10) 

uo > a in Qq, uq < a in Qq, (l-H) 

where denotes the region enclosed by Fq and the region enclosed 
between and Fq. 

It is standard that problem (P^) has a unique smooth solution, which we 
denote hy u^. As e ^ 0, a formal asymptotic analysis shows the following: 
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in the very early stage, the diffusion term An is neghgible compared with 
the reaction term e~'^{f{u) — eg"^ {x,t,u)) so that, in the rescaled time scale 
r = t/e'^, the equation is well approximated by the ordinary differential 
equation Ur = f{u) + 0{e). Hence, in view of the profile of /, the value of tt^ 
quickly becomes close to either a_|_ or a„ in most part of fi, creating a steep 
interface (transition layer) between the regions {u^ ~ a_} and {u^ ~ Oi+}- 
Once such an interface develops, the diffusion term becomes large near the 
interface, and comes to balance with the reaction term. As a result, the 
interface ceases rapid development and starts to propagate in a much slower 
time scale. 

To study such interfacial behavior, it is useful to consider a formal asymp- 
totic limit of (P^) as e — > 0. Then the limit solution u{x,t) will be a step 
function taking the value a+ on one side of the interface, and a_ on the 
other side. This sharp interface, which we will denote by Tt, obeys a certain 
law of motion, which is expressed as follows (see Section [2] for details): 



\t=o 



-{N — 1)k + co{G{x, t, «+) — G{x, t, «-)) on Ft 

= ro. 



where Vn is the normal velocity of Tt in the exterior direction, k the mean 
curvature at each point of Ft, 



Co 



fO'+ -1—1 

V2 {W{s) -W{a-))^/^ds , (1.12) 

J a- 



W{s) = - f{r)dr, Gix,t,s)= g{x,t,r)dr. 

J a J a 

It is well known that problem {P^) possesses locally in time a unique smooth 
solution. Let < t < T^°-^ , j<max g (Q, +oo], be the maximal time interval 
for the existence of the solution of {P^) and denote this solution by F = 
Uo<f<T'""^ (-'^i X {i}). Hereafter, we fix T such that < T < T*""^ and work 
on [0,T]. More precisely, so as g{-,-,u), the function G(-,-,n) is of class 
^1+1?, -2- ^ which implies, by the standard theory of parabolic equations, 
that r is of class C^^"^—. For more details, we refer to [9], Lemma 2.1. 
Next we set 

Qt:=^x (0,r), 

and, for each t G [0,T], we denote by the region enclosed by the hyper- 
surface Tt, and by Qf' the region enclosed between dQ and Ft. We define a 
step function u{x,t) by 

I a+ in nf , , , , 

u{x,t) = { + . * fortG[0,r], (1.13) 
I a_ m 
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which represents the formal asymptotic hmit of (or the sharp interface 
limit) as e — > 0. 

The aim of the present paper is to make a detailed study of the limiting 
behavior of the solution of problem (P^) as e — > 0. Our first main result, 
Theorem 11.31 describes the profile of the solution after a very short initial 
period. It asserts that: given a virtually arbitrary initial data uq, the solu- 
tion quickly becomes close to a±, except in a small neighborhood of the 
initial interface Tq, creating a steep transition layer around Tq {generation 
of interface). The time needed to develop such a transition layer, which 
we will denote by t^ , is of order lne|. The theorem then states that the 
solution u'^ remains close to the step function u on the time interval [i^,?^] 
(motion of interface); in other words, the motion of the transition layer is 
well approximated by the limit interface equation (P^)- 

Theorem 1.3 (Generation and motion of interface). Let r] be an arbitrary 
constant satisfying < rj < min(a — a_ , a+ — a) and set 

fi = f'{a). 

Then there exist positive constants Eq and C such that, for all e G (0,eo) 
and for all t^ < t <T, where t^ := fi~^e'^\ lne\, we have 

[a_ -r],a+ + 7]] if x £ MceO^t) 
u%x, t) G - 77, a_ + 7?] if x € \ MceiTt) (1.14) 

, [a+ -r],a+ + ri] if x € \ MceO^t), 



where J\fr{Tt) '■= {x G il.,dist{x,Tt) < r} denotes the r -neighborhood ofT- 



Corollary 1.4 (Convergence). As e ^ 0, converges to u everywhere in 

The next theorem is concerned with the relation between the actual 
interface Ff := {x G O, u^{x,t) = a} and the formal asymptotic limit Ft, 
which is given as the solution of (P^). 

Theorem 1.5 (Error estimate). There exists C > such that 

rf C Mceirt) for 0<t<T. (1.15) 

Corollary 1.6 (Convergence of interface). There exists C > such that 

dnm,'^t) <Ce for 0<t<T, (1.16) 

where d-}i{A,B) := maxjsup^g^ (i(a, S), supfj^^ d{b, A)} denotes the Haus- 
dorff distance between two compact sets A and B. Consequently, Ff — > Vt 
as e ^ uniformly in < t <T, in the sense of Hausdorff distance. 
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Note that the estimates (jl.lSp and ()1.16p follow from Theorem 11.31 in 
the range < t < T, but the range < t < has to be treated by a 
separate argument since the behavior of the solution in this time range is 
quite different from that of the later stage. 

The estimate (I1.14p in our Theorem 11.31 implies that, once a transition 
layer is formed, its thickness remains within order e for the rest of time. 
Here, by "thickness of interface" we mean the smallest r > satisfying 

{ X e 0, u'^(j;, t) [a- - 7], a- +T]]U [a+ - t], a+ + r]]} C Afr{Tl). 

Naturally this quantity depends on rj, but the estimates ()1.14p and ()1.16p 
assert that it is bounded by 2Ce (with the constant C depending on i]) 
regardless of the choice of r/ > 0. 

Remark 1.7 (Optimality of the thickness estimate). The above 0{e) estimate 
is optimal, i.e. the interface cannot be thinner than this order. In fact, 
rescaling time and space as t := t/e'^, y := x/e, we get 

Ur = AyU + f{u) - e/. 

Thus, by the uniform boundedness of u and by standard parabolic estimates, 
we have |Vj^u| < M for some constant M > 0, which implies 

\^xu{x,t)\ < —. 

From this bound it is clear that the thickness of interface cannot be smaller 
than M~^{a^ — a_) e, hence, by (jl.l4p . it has to be exactly of order e. Intu- 
itively, this 0(e) estimate follows also from the formal asymptotic expansion 
(j2.3p , but the validity of such an expansion is far from obvious for solutions 
with arbitrary initial data. □ 

Our 0(e) estimate is new, even in the special case where = 0, pro- 
vided that N > 2. Previously, the best thickness estimate in the literature 
was of order e|lne| (see [7j), except that X. Chen has recently obtained 
an 0(e) estimate for the case N = 1 by a different argument (private 
communication). We also refer to the forthcoming papers [20] and |19j . 
in which the same 0(e) estimate is established for different but related 
problems. The paper [20] is concerned with a "balanced type" Allen-Cahn 
equation with large spatial inhomogeneity, namely an equation of the form 
ut = V{k{x)\/u) + e~'^ h{x) f (u) , and [19] is concerned with a Lotka-Volterra 
competition-diffusion system with large spatial inhomogeneity whose non- 
linearity is of the balanced bistable type. 

Remark 1.8 (Optimality of the generation time). The estimate ()1.14p also 
implies that the generation of interface takes place within the time span 
of t^. This estimate is optimal. In other words, a well-developed interface 
cannot appear much earlier; see Proposition 13.101 for details. □ 
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The singular limit of the Allen-Cahn equation was first studied in the 
pioneering work of Allen and Cahn [2] and, slightly later, in Kawasaki and 
Ohta [23] from the point of view of physicists. They derived the interface 
equation by formal asymptotic analysis, thereby revealing that the interface 
moves by the mean curvature. These early observations triggered a flow of 
mathematical studies aiming at rigorous justification of the above limiting 
procedure; see, for example, [B], [ZIE] and [MIES] for results in the frame- 
work of classical solutions, and [E], [H [5] and [21] for the case where Ff is 
a viscosity solution of the interface equation. 

As for problem (P^), whose nonlinearity is slightly unbalanced, the limit 
interface equation involves a pressure term as well as the curvature term as 
indicated in (-P"). This fact has been long known on a formal level; see e.g. 
|27| . Ei, lida and Yanagida [13] proved rigorously that the motion of the 
layers of (.P^) is well approximated by the limit interface equation {P^), on 
the condition that the initial data has already a well developed transition 
layer. In other words, they studied the motion of interface, but not the 
generation of interface. 

1.2 Singular limit of reaction-diffusion systems 

Our results can be extended to reaction-diffusion systems of the form 

in O X (0, -l-oo) 
in J7 X (0, +oo) 

on 9r2 X (0, +oo) 

in Q 
in Q, 

where D is a positive constant, and f^, h are functions such that 

(F) there exist functions fi{u,v), f2{u,v) such that 

fiu,v) = fiu) + ehiu,v) + e^fi{u,v), (1.17) 

where f{u) is a bistable nonlinearity satisfying (jl.l|) . (|1.2|) . and /I, 
along with its derivatives in u, v, remain bounded as e ^ 0; 

(H) for any constant L,M > there exists a constant Mi > M such that 

h{u,-Mi) >0> h{u,Mi) for \u\ < L. (1.18) 

The conditions (F) and (H) imply that the ODE system 



ut = Au + ^ f{u,v) 

vt = DAv + h{u, v) 
du dv ^ 
du du 
u{x, 0) = uo{x) 

v{x,0) = vo{x) 



u=—^f^{u,v), v = h{u,v) 
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has a family of invariant rectangles of the form {\u\ < L, \v\ < M}, provided 
that e is sufficiently small. The maximum principle and standard parabolic 
estimates then guarantee that every solution of (RD^) exists globally for 
t > and remains bounded as t ^ oo (see Proposition I7.ip . Apart from 
(jl.lSp . we do not make any specific assumptions on the function h. 

Problem {RD^) represents a large class of important reaction-diffusion 
systems including the FitzHugh-Nagumo system 



^u + ^{f{u)-ev) 



vt = DAv + au — Pv, 



(1.19) 



which is a simplified model for nervous transmission, and the following type 
of prey-predator system: 



Ut = Au + ^{{1 -u){u- 1/2) -ev)t 
Vt = DAv + {au — j3v)v. 



(1.20) 



Remark 1.9. In some equations such as the prey-predator system (ll.20p . 
only nonnegative solutions are to be considered. In such a case, we replace 
the condition (|1.18p by 

h{u, 0) > > /i(n, Ml) for < n < L, 

and assume /^(O, v) > 0. The rest of the argument remains the same. □ 

Now the same formal analysis as is used to derive (-P") shows that the 
singular limit of {RD^) as e — > is the following moving boundary problem: 



{RD^) 



Vn = -{N -l)K-CQFi{v{x,t)) 

Vt = DAv + h{u, v) 



^4=0-^0 



1^ = 



v{x, 0) = vq{x) 
where u is the step function defined in (I1.13P and 

= / fi{r,v) dr. 



on Tt 

in O X (0, +oo) 

on di} X (0, +00) 

in 0,, 



This is a system consisting of an equation of surface motion and a partial 
differential equation. Since u is determined straightforwardly from Tt, in 
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what follows, by a solution of (RD^) we mean the pair (F, v) := {Tt, v{x, t)). 
In the case of the FitzHugh-Nagumo system ()1.19p . (RD^) reduces to 

Vn = —{N — 1)k + Co (a+ — a-)v{x, t) 

Vt = DAv + au — (3v, 

while in the prey-predator system (ll.20p . (RD^) reduces to 

Vn = -{N -l)K + coi{x,t)/2 

Vt = DAv + (an — (3v)u. 

Note that the positive sign in front of the term CQv{x,t) in the interface 
equation implies an inhibitory effect on u, since the velocity Vn is measured 
in the exterior normal direction, toward which u decreases. 

Lemma 1.10 (Local existence). Assume that vq £ (7^(0) and that Tq is a 
hypersurface which is the boundary of a domain Dq CC O. Then there 
exists T"*"^' S (0, +00] such that the limit free boundary problem (RD^) has 
a unique solution {T,v) in the interval [0, T™'"^). 

This existence result was established in [10]. The uniqueness can be 
obtained by using the estimates in [S] . 

Hereafter, we fix T such that < T < T"*''^ and work on [0,r]. Our 
main results for the system (RD^) are the following: 

Theorem 1.11 (Thickness of interface) . Let (11.170 and (jl.lSp hold (or let 
the assumptions in Remark \l.y\ hold). Assume also that uq satisfies (jl.lUp 
and (jl.lip . Then the same conclusion as in Theorem \1.3\ holds for {RD'^). 

Corollary 1.12 (Convergence). Under the assumptions of Theorem \1.11\ 
the same conclusion as in Corollary \1.4\ holds for {RD^). 

Theorem 1.13 (Error estimate). Let the assumptions of Theorem \ 1 . ll\ hold. 
Then the same conclusion as in Theorem li.5l holds for [RD^). Moreover, 
there exists a constant C > such that 

ll'ii^ - 'i'llL°°(nx(o,T)) < Ce. 

Corollary 1.14 (Convergence of interface). Under the assumptions of The- 
orem [Tm[ the same conclusion as in Corollary \1.6\ holds for {RD^). 

The organization of this paper is as follows. In Section [21 we derive the 
interface equation {P^) from (-P^) by formal asymptotic expansions which 
involve the so-called signed distance function. 

In Sections [3| and [U we present basic estimates concerning the generation 
of interface for {P^). For the clarity of underlying ideas, we first consider 
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the special case where 5^ = in Section [Sj and deal with the general case in 
Section HI 

In Section [5] we prove a preliminary result on the motion of interface 
(Lemma I5.ip . which implies that if the initial data has already a well- 
developed transition layer, then the layer remains to exist for < i < T 
and its motion is well approximated by the interface equation (P^)- 

Our approach in Sections O to [S] is based on the sub- and super-solution 
method, but we use two completely different sets of sub- and super-solutions. 
More precisely, the sub- and super-solutions for the motion of interface are 
constructed by using the first two terms of the formal asymptotic expansion 
(j2.3p , while those for the generation of interface are constructed by modifying 
the solution of the equation in the absence of diffusion: ut = e~^/(u). 

In Section [6l we prove our main results for (P^): Theorems ll.3| [T3] and 
their respective corollaries. 

In the final section, we study the reaction-diffusion system (RD^) and 
prove Theorems 11.111 fl.l3l and their corollaries. These results are obtained 
by applying a slightly modified version of the results for (P^)- The strategy 
is to regard f^{u,v) as a perturbation of f{u). Indeed, the equation for u 
in (RD^) is identical to (P^) if we set = — /i — However, what 

makes the analysis difficult is the fact that is no longer a given function 
but a quantity that depends on the unknown function v'^. In particular, 
the existence of the limit g^ ^ g {e ^ 0) is not a priori guaranteed, and 
the estimate ()1.6p is far from obvious. As it turns out, the standard or 
Schauder estimates for would not yield (|1.6p . because of the fact that 
converges to a discontinuous function as e — > 0. In order to overcome this 
difficulty, we derive a fine estimate of that is based on estimates of the 
heat kernel and the fact that remains uniformly smooth outside of an 
0(e) neighborhood of the smooth hypersurface Fj. 

2 Formal derivation of the interface motion equa- 
tion 

In this section we derive the equation of interface motion corresponding to 
problem (P^) by using a formal asymptotic expansion. The resulting inter- 
face equation can be regarded as the singular limit of (P"^) as e — > 0. Our 
argument is basically along the same lines with the formal derivation given 
by Nakamura, Matano, Hilhorst and Schatzle [26], who studied a similar 
but slightly different type of spatially inhomogeneous equations by formal 
analysis. Let us also mention some earlier papers [1], [T7j and [27] involving 
the method of matched asymptotic expansions for problems that are related 
to ours. 

As in [26], the first two terms of the asymptotic expansion determine 
the interface equation. Though our analysis in this section is for the most 



10 



part formal, the observations we make here wih help the rigorous analysis 
in later sections. 



Let be the solution of (-P^). We recall that Ff := {x G ^l, u^{x, t) = a} 
is the interface at time t and call := IJf>o(-'^t ^ i^i) interface. Let 
r = \jQ<t<T^t ^ {^}) unique solution of the limit geometric motion 

problem {P^) and let d be the signed distance function to T defined by: 

I — dist(x,rt) for X G , 

where dist(x, F^) is the distance from x to the hypersurface Ft in 0. We 
remark that d = on F and that |Vd| = 1 in a neighborhood of F. We then 
define 

0<t<T 0<t<T 

We also assume that the solution is of the form 

u%x,t) = a± + euf{x,t) -\ in (2.2) 

away from the interface F (the outer expansion), and 

u'{x,t) = Uo{x,t,0 + eUi{x,t,0 + ■ ■ ■ (2.3) 

near F (the inner expansion), where Uj{x, t, z), j = 0, 1, • • • , are defined for 
xGfi, t>0, zGM and ^ := d{x,t)/e. The stretched space variable ^ gives 
exactly the right spatial scaling to describe the rapid transition between the 
regions {ti*" ~ a_} and {u^ w Q;+}. We normalize Uq in such a way that 

Uo{x, t,0) = a 

(normalization conditions). To make the inner and outer expansions consis- 
tent, we require that 

Uo{x,t,+oo) = a+, Uo{x,t,—oo) = a-. (2-4) 

As we will see below this will determine Uq uniquely, which will then deter- 
mine Ui. 

In what follows we will substitute the inner expansion (|2.3|) into the 
parabolic equation of problem (-P^) and collect the and terms. To 
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that purpose we compute the needed terms and get 

ul = Uot + Uozj + eUu + UiJt + ■■■ 
Vd 

Vn^ = VUq + Uoz h eVUi + U^Vd H 

e 

Au' = AUo + 2— • Vf/o. + Uoz— + Uozz^-^^ + eAUi 
e e e'^ 

iVcii^ 

+2Vd • Vf/i^ + UizAd + C/i^^J ^ + • • • 

e 

/(n") =/(f/o)+£/'(t/o)f/i + 0(£') 

/(x,*,^^) = c/(2;,t,ii^) +0(e) ( < — in view of ([rn]) ) 
= 5(x,t,[/o) +0(e), 

where the functions Ui [i = 0, 1), as well as their derivatives, are taken at 
the point {x,t,d{x,t)/e). Note also that V and A stand for V^; and Ax, 
respectively. Collecting the terms yields 

Uozz + f{Uo) =0. 

In view of the normalization and matching conditions, we can now assert 
that Uo{x, t, z) = Uo{z), where Uo{z) is the unique solution of the stationary 
problem 

f uo" + fm = 

I C/o(-oo) = a_, [/q(0) = a, Uq{+oo) = a+. 

This solution represents the first approximation of the profile of a transition 
layer around the interface observed in the stretched coordinates. Note that 
the integral condition ()1.2p guarantees the existence of a solution of (12. 5p . 
For example, in the simple case where f{u) = u{l — u^), we have Uo{z) = 
tanh(^;/\/2). In the general case, the following standard estimates hold: 

Lemma 2.1. There exist positive constants C and A such that 

< a+ - Uq[z) < Ce-^l^l forz>0 
< Uo{z) - a_ < Ce-^l^l for z<0. 

In addition, Uq is a strictly increasing function and, for j = 1,2, 

\DWoiz)\ < Ce~^l^l for z G M. (2.6) 
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Proof. We only give an outline. Rewriting the equation in (12. 5p as 

u = v, v = -f{u), 

we see that {Uo{z),Uq{z)) is a heteroclinic orbit of the above system con- 
necting the equilibria (q_, 0) and (a+, 0). These equilibria are saddle points, 
with the linearized eigenvalues {A_, — A_} and {A+, — A+j, respectively, 
where 

A_ = ^-/'(a_), A+ = 
Consequently, we have 

, , _ J a_ +Cie^-^ + o(e^-^) as z ^ -oo, , . 

~ \ Q+ + C2e-^+^ + o(e-^+^) asz^+oo, ^ ' 

for some constants Ci , C2 . The desired estimates now follow by setting 
A = min(A+, A_). □ 

Next we collect the terms. Recalling that VUqz = and that |V(i| = 
1 near Tf, we get 

Ui,, + f'mUi = Uo'idt - Ad) + g{x, t, Uo). (2.8) 

This equation can be seen as a linearized problem for (j2.5p with an inho- 
mogeneous term. As is well-known (see, for instance, |26]). the solvability 
condition for the above equation plays the key role in determining the equa- 
tion of interface motion. The following lemma is rather standard, but we 
give an outline of the proof for the convenience of the reader. 

Lemma 2.2 (Solvability condition). Let A{z) be a bounded function on 
—00 < z < 00. Then the problem 



( ^P,, + f'{Uo{z))iP = A{z), ze 
\ Tp{0) =0, ^£ L°°(]R), 
has a solution if and only if 



(2.9) 



A{z)Uo\z)dz = 0. (2.10) 

Moreover the solution, if it exists, is unique and satisfies 

|V'(z)| < CPIIioo /orzGM, (2.11) 
for some constant C > 0. 
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Proof. Multiplying the equation by Uq and integrating it by parts, we 
easily see that the condition (j2.10p is necessary. Conversely, suppose that 
this condition is satisfied. Then, since Uq is a bounded positive solution to 
the homogeneous equation il)zz + f iUo{z))ij; = 0, one can use the method of 
variation of constants to find the above solution ijj explicitly. More precisely. 



Jo J^oo .2.12) 

i-z PCX) y ^ 



where ip := Uq . The estimate ()2.1ip now follows from the above expression 
and (EITD. □ 



dz = 0, 



From the above lemma, the solvability condition for (12. 8p is given by 
Uo'\z){dt-^d){x,t) + g{x,tMz))Uo'{z) 
for all G Qt- Hence we get 

d - Ad= /Rg(^'^'^o(^))^o'(^) '^^ 
J^Uo'^{z)dz 

which gives 

~ ^ G{x, t,a+) - G{x, t, a-) 



J^Uo'^iz) dz 

Moreover, multiplying equation (j2.5p by Uq' and integrating it from —00 to 
z, we obtain 

= r {Uo"Uo' + f{Uo)Uo'){s)ds 



= \uo'\z)-W{Uo{z)) + W{a^), 

where we have also used the fact that C/o(— 00) = q_ and C/o'(— 00) = 0. 
This implies that 

Uo'{z) = V2{WiUo{z))-W{a.)f\ 

and therefore 

Uo''^iz)dz = [ Uo'{z)V2{WiUoiz)) -W{a.)f^^dz 
Jr 

= V2 [ ''{W{s) - W{a-))^^^ds. 
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It then follows, in view of the definition of cq in ()1.12p . that 

dt = Ad- co{G{x,t,a+) - G{x,t,a-)). (2.13) 

We are now ready to derive the equation of interface motion. Since Vd (= 
V^; d{x, t)) coincides with the outward normal unit vector to the hypersurface 
Tt, we have dt{x,t) = —Vn, where Vn is the normal velocity of the interface 
r^. It is also known that the mean curvature k of the interface is equal to 
Ad/{N — 1). Thus the equation of interface motion is given by: 

Vn = -iN -l)K + coiG{x,t,a+)-G{x,t,a^)) on Tt. (2.14) 

Summarizing, under the assumption that the solution of problem (P^) 
satisfies 

I a+ in 

u ^ < _ as e ^ 0, 

I a_ in Qrp, 

we have formally proved that the boundary Tt between and moves 
according to the law (j2.14p . 

To conclude this section, we give basic estimates for Ui{x, t, z), which we 
will need in Section \5\ to study the motion of interface. Substituting ()2.13p 
into (|2.8p gives 



Ui,, + f'{Uo{z))Ui = gix,t,Uo{z)) - j{x,t)Uo'{z) 
Ui{x,t,0)=0, Ui{x,t,-)eL' 



(2.15) 



where 

7(x,t) = co{G{x,t,a+) - G{x,t,a-)). (2.16) 
Thus Ui{x,t,z) is a solution of ()2.9p with 

A = Ao{x, t, z) := g{x, t, Uo{z)) - 7(x, t)Uo'{z), (2.17) 

where the variables x, t are considered parameters. The problem ()2.15p has a 
unique solution by virtue of Lemma [2.21 Moreover, since Ao{x, t, z) remains 
bounded as (x, t, z) varies in $7 x [0, T] x M, the estimate ()2.1ip implies 

\Ui{x,t,z)\ <M for X G t G [0,r], z G M, (2.18) 

for some constant M > 0. Similarly, since Vt/i is a solution of ()2.9p with 

A = \I^Ao{x,t,z) (^ = V4g{x,t,Uo{z)) --f{x,t)Uo'{z)) ), 

and since g is assumed to be G^ in x, we obtain 

\V^Ui{x,t,z)\ < M for X GO, t G [0,r], z G M, (2.19) 
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for some constant M > 0. 

To obtain estimates as z — > ±00, we first observe that (j2.7p implies 

Ao{x,t,z) - g{x,t,a±) = 0{e~^\''\) as z ^ ±00, (2.20) 

uniformly in x S 0, t G [0,7"]. We then apply the following general esti- 
mates: 

Lemma 2.3. Let the assumptions o f Lemma \2. S\ hold, and assume further 
that A{z) — = 0(e~^l^l) as z ^ ±00 for some constants A'^ and 
(5 > 0. Then there exists a constant A > such that 

^i') - = 0{e~'\^\), \^p'{z)\ + \^p"{z)\ = 0{e~'\^\), (2.21) 

/ (a±) 

as z — > ±00. 

Proof. We only state the outline. To derive the former estimate, we need 
a slightly more elaborate version of ()2.7p . Since f{u) is C^, we have f{u) = 
{u — Oi±)f'{a±) + 0((u — a-t)^) . Consequently, 

_ r a_ + Ci e^-" + 0(e2^-^) as z ^ -00, 

a+ + C2e-^+^ + 0(e-2^+-) as 2 ^ +00. ^^'^^^ 



Using the expression (j2.12p along with the estimate A{z) — A^ = 0{e "^'^1) 
and (I2.22p . we see that 

_A^ 

(Ad 



^p(z) = -7^^ + 0{ |z|e-^±l"l ) + o(e-™'^('5'^±)l^l ) as z ^ ±cx). 



This implies the former estimate in (I2.2ip . where A can be any constant sat- 
isfying < A < min(A_, A+, 5). Substituting this into equation (12. 9p gives 
the estimate for ipzz- Finally, the estimate for ipz follows by integrating il^^z 
from ±00 to z. □ 

Prom the above lemma and ()2.20p we obtain the estimate 

\Uu{x,t,z)\ + \Ui,,ix,t,z)\ <Ce^^l^l, (2.23) 

for x € fl, t £ [0, T], 2: G ]R. Similarly, since ()2.6p implies 

{V^Ao){x,t,z) - {V^g){x,t,a±) = 0{e-^\'\) as z ^ ±00, 

we can apply Lemma 12.31 to -0 = V^f^i, to obtain 

\V^^Ui,{x,t,z)\ + \V^Uu,{x,t,z)\ < C7e-^l^l 

for X E 0, t G [0, T], z S M. As a consequence, there is a constant, which we 
denote again by M, such that 

\V,Uu{x,t,z)\<M. (2.24) 
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Next we consider the boundary condition. Note that ()1.7p implies 
—Ao = -^[g{x,tMz))-^{x,t)Uo'{z)\=Q on OQ. (2.25) 
Consequently, from the expression (|2.12|) . or equivalently the expression 

Ui (x, t, z) = U'^{z) ( (f/^ (0) y (x, t, C/o (0 di) dC, 
we see that 

— ^ = on dn. (2.26) 

3 Generation of interface: the case ^^^^ = 

This section deals with the generation of interface, namely the rapid for- 
mation of internal layers that takes place in a neighborhood of Tq = {x £ 
Q, uo{x) = a} within the time span of order lne|. For the time being we 
focus on the special case where = 0. We will discuss the general case in 
Section HI In the sequel, rjo will stand for the following quantity: 

r/o := min(a — a_, a+ — a). 

Our main result in this section is the following: 

Theorem 3.1. Let rj G (OjTJq) be arbitrary and define fi as the derivative of 
f{u) at the unstable zero u = a, that is 

fx = f'{a). (3.1) 

Then there exist positive constants Eq and Mq such that, for all e £ (0,eo), 

(i) for all X £ Q, 

a- — rj < u'^{x, fJ,~'^e'^\ lne|) < a+ + rj, (3.2) 

(ii) for all x £ 0, such that |no(x) — a\ > Mqe, we have that 

if uo{x) > a + Mqe then u'^ {x, fi~^e'^\lne\) > a+ — r], (3.3) 
if uo{x) < a — Mqe then u"^ {x, fi~^E'^\lnE\) < + r]. (3.4) 

The above theorem will be proved by constructing a suitable pair of sub- 
and super-solutions. Note that we do not need condition (11.21) in proving 
this theorem. 
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3.1 The bistable ordinary differential equation 

Let us first consider the problem without diffusion: 

1 



ut 



f{u), u{x,0) = uo{x). 



This solution is written in the form 



u{x,t) = y(^-^, uo{x)J, 
where Y{t, ^) denotes the solution of the ordinary differential equation 



YriT,0 =/(y(r,0) forr>0 

me) =e 



(3.5) 



Here ^ ranges over the interval (— 2Co,2Co), with Co being the constant 
defined in (jl.Sp . We first study basic properties of Y. 



Lemma 3.2. We have > 0, for all ^ ^ {a-, a, a+}, r > 0. Furthermore, 



Ydr,0 



Proof. First, differentiating equation (|3.5p by ^, we obtain 

f Y^, = Yl:f{Y) 

1 n(o,e) = i, 

which is integrated as follows: 

Y^{t,0=^^v \ r f{Y{s,i))ds >0. 



(3.6) 



We then differentiate equation ()3.5p by r and obtain 

f Yrr = Yrf'iY) 
\ Yr{0,0=f{0, 

which in turn implies 

YriT,0 =f{Oew\rf'iYis,0)ds 



= f{OYdr,C)- 

This last equality, in view of (j3.5p , completes the proof of Lemma 13.21 □ 
For ^ ^ {«_, a, a+}, we define a function A(t, ^) by 

f'(Y{T,o)-no 



(3.7) 
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Lemma 3.3. We have, for all ^ ^ {«-, a, a+}, r > 0, 

A{t,C)= r f"{Y{s,C))Y^{s,^)ds. 
Jo 

Proof. Differentiating by ^ the equality of Lemma 13.21 leads to 

% = yl(r,e)n, (3.8) 
whereas differentiating (13. 6p by ^ yields 

Jo 

These two last results complete the proof of Lemma 13.31 □ 

Next we need some estimates on Y and its derivatives. First, we estimate 
the speed of the evolution of Y when the initial value ^ lies between a_ + r] 
and Q+ — Tj. 

Lemma 3.4. Letr] G (0,%) be arbitrary. Then there exist positive constants 
Ci = Ciiji), C2 = 02(11) <^'^d C3 = C-iij]) such that 

(i) if ^ ^ (a, a+ — rj) then, for every r > such that Y{t, remains in 
the interval (a, a;+ — r]), we have 

Cie^'^ <Yi:{T,0<C2e>'\ (3-9) 

|^(T,e)| <C3(e^--l), (3.10) 
where fi is the constant defined in (|3.ip .' 

(ii) if ^ & (a- +r],a) then, for every r > such that Y{t,(,) remains in 
the interval (a_ + 7/, a), we have (13. 9p and ()3.10p . 

Proof. We take ^ G (a,a+ — rj) and suppose that, for s £ (0,r), Y{s,S,) 
remains in the interval (a, — r/). Integrating the equality 



fiYis,0) 

from to r yields 



lo fiyis,o) 

Hence by the change of variable q = Y{s,(^) we get 
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Moreover, the equality of Lemma 13.21 leads to 

-y(r,o f>i 



^(-'^/(a) , ni)-n^) u, (3 13) 



/^r + y h{q)dq 



where h{q) = {f'{q) — fi)/f{q). As h{q) tends to f"{a)/ f'{a) when q tends 
to a, h is continuous on [a, a+ — r/] . Hence we can define 

H = H{7]) := \\h\\L^(^a,a+-ri)- 

Since |y(r, ^) — takes its value in the interval [0, «+ — a — tj] C [0, a+ — a], 
it follows from ()3.13p that 

fiT — H{a+ — o) < In l^(r, ^) < fiT + H{a^ — a), 

which, in turn, proves ()3.9p . Next Lemma 13.31 and ()3.9p yield 

JO 

< C3(e^- - 1), 

which completes the proof of ()3.10p . The case where ^ and y(r, ^) are in 
(a_ + rj, a) is similar and omitted. □ 

Corollary 3.5. Let r] G (0,170) be arbitrary. Then there exist positive con- 
stants Ci = Ci(r/) and C2 = C2{r]) such that 

(i) if C ^ (a,a+ — rj) then, for every r > such that y(r, ^) remains in 
the interval (a, q;+ — rj), we have 

Cie'^"(e - a) < y(r,0 - a < C72e^"(e - a); (3.14) 

(ii) if i ^ (o- + a) t/ien, /or every r > suc/i i/iai ^(r, 1^) remains in 
the interval (a_ +ry, a), w;e /laue 

C2e'^^(,t - a) < y(r,0 - a < Cie^^(^ - a). (3.15) 

Proof. We can find Bi = Bi{rj) > and B2 = B2{r]) > such that, for all 
q G (a,a+ - ??), 

Bi{q-a)< f{q)<B2{q-a). (3.16) 
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We use this inequality for a < y(r, < a+ — t] to obtain 

Bi{Y{t,0 - a) < f{Y{T,0) < B2{Y{t,0 - a). 
We also use this inequality for a < ^ < a+ — to obtain 

Bi{^-a)<f{0<B2{^-a). 
Next we use the equality 1^ = f{Y)/f{^) of Lemma 13.21 to deduce that 

|i(y(r,0 - a) < (e - a)Y^{r,0 < ^{Y{t,0 " «), 



which, in view of ()3.9p . implies that 

fiCie^-le - a) < y(r,e) - a < ^C2e^^iC - a). 

02 -Dl 

This proves (|3.14|) . The proof of (j3.15p is similar and is omitted. □ 

We now present estimates in the case where the initial value ^ is smaller 
than a_ + r/ or larger than a_|_ — rj. 

Lemma 3.6. Let rj G (0, ryo) and M > be arbitrary. Then there exists a 
positive constant C4 = C^^rjjM) such that 

(i) if i ^ [a+ — a+ + M], then, for all r > 0, Y{t,£,) remains in the 
interval [a-|_ — r/, a+ + M] and 

|A(r,OI <C4r for r>0; (3.17) 

(ii) if ^ ^ [«- — M, a_ + rj], then, for all r > 0, Y{t,S,) remains in the 
interval [a_ — M, a- + r]] and (j3.17p holds. 

Proof. Since statement (i) and statement (ii) can be treated in the same 
way, we will only prove the former. The fact that Y{t, remains in the 
interval [a+ — ry, a+ + M] directly follows from the bistable properties of /, 
or, more precisely, from the sign conditions /(a+ — ry) > 0, /(a+ + M) < 0. 

To prove ()3.17p . suppose first that ^ G [a+,a+ + M]. In view of (jl.ip . /' 
is strictly negative in an interval of the form [a^ , a_|_ + c] and / is negative in 
[q:+, 00). We denote by —m < the maximum of / on [a+ + c, M]. Then, as 
long as Y{t, ^) remains in the interval [a+ + c, M], the ordinary differential 
equation (13. 5p implies 

Yr < —m. 

This means that, for any G [a+, a+ + M], we have 

M-c 



Y{t,S^) G [a+,a+ + c] for t > t :- 



m 
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In view of this, and considering that f'{Y) < for Y G + c], we see 

from the expression (j3.6p that 



Yd^,^) =exp 

< exp 

< exp 



f'{Yis,0)ds 
nYis,C))ds 



exp 



f'iY{s,0)ds 







sup \f'{z)\ds 

^JO zG[o--M,a++M] 



= : C4, 



for all T > r. It is clear from the same estimate p. 60 that Y^ < C4 holds 
also for < r < r. We can then use Lemma 13.31 to deduce that 



|^(r,OI <^4 / \riY{s,mds<C^T. 



The case ^ € [a+ — t], a+] can be treated in the same way. This completes 
the proof of the lemma. □ 

Now we choose the constant M in the above lemma sufficiently large so 
that [-2Co, 2Co] C [a_ - M,a+ + M], and fix M hereafter. Then C4 only 
depends on rj. Using the fact that r = 0(e^^ — 1) for r > 0, one can easily 
deduce from ()3.10p and (13.170 the following general estimate. 

Lemma 3.7. Let r] £ (0, 770) be arbitrary and let Cq be the constant defined 
in (jl.Sp . Then there exists a positive constant C5 = C^rj) such that, for all 
T>0 and all ^ G (-2Co, 2Co), 



|^(r,e)| <C5(e'^^-l) 



3.2 Construction of sub- and super-solutions 

We are now ready to construct the sub- and super-solutions for the study 
of generation of interface. For simplicity, we first consider the case where 

^ = on dn. (3.18) 

In this case, our sub- and super-solutions are given by 

wf{x,t) = ^(^, uo{x)±e'Ce{e^'/'" - 1)). (3.19) 

In the general case where (13.180 does not necessarily hold, we have to slightly 
modify wf{x,t) near the boundary dQ. This will be discussed later. 
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Lemma 3.8. Assume (jS.lSp . Then there exist positive constants Eq and 
Cq such that, for all e G (0, eo)? {''^e^'^t) "^^ ^ P^''''^ ^'^'^ super- 

solutions for problem (P^), in the domain Q x [0, /i~"'^e^| lne|], satisfying 
w~{x,0) = w'^{x,0) = uq{x). Consequently 



{x,t) <u^{x,t) <wf{x,t) for X £ <t < 1^1 ^e^|lne| 
Proof. The assumption (IS.lSp implies 

-^ = on X (0,+oo). 
Now we define an operator Cq by 



Cqu := ut- Au ^f{u), 

and prove that Cqw^ > . Straightforward computations yield 



(3.20) 



' £2- 



therefore, in view of the ordinary differential equation (|3.5p . 



Cowl. 



MA^_Ano-^|Vno| 



Y,. 



We note that, in the range < t < fi ^e^j lne\, we have, for eq sufficiently 
small, 

< e^Ceie^'/'" - 1) < e^Ceie^' - 1) < Co, 



where Co is the constant defined in (II. Sp . Hence 

e := uo{x) ± C^{e^'l'" - 1) G (-2Co, 2Co), 

and it follows from the estimate of ^4 = Y^^jY^ in Lemma 13.71 with the 
choice r := tje^ , that 



> 



Cefie^'/'" - lAnol - C^ie^'/'" - 1)\Vuo\^\y^ 
"(Ce /i - Csl Vtxo|2)e^*/^' - |Ano| + C5I VuqP 
Since 1^ > 0, this inequality implies that, for Cq large enough, 



> 



Cowt > 



Cafi — C^Co^ — Co 



Y^ > 0. 



Hence is a super-solution for problem (P^). Similarly it;^ is a sub- 
solution. Obviously 0) = w'^{x,0) =uq{x). Lemma[3]8]is proved. □ 
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In the more general case where (I3.18P is not necessarily valid, one can 
proceed as follows: in view of (jl.lOp and (jl.lip there exist positive constants 
di, p such that uo{x) > a + p if d{x,dO,) < di. Let x be a smooth cut-off 
function defined on [0, +00) such that < x 1) x(0) = x'(0) = and 
x{z) = 1 for z > di. Then we define 

Uq{x) = x{d{x, d^)) uq{x) + [1 — x{d{x, max no(x), 

Uq{x) = x{d{x,dn))uo{x) + [1 - x{d{x,dfl))] (a + p). 

Clearly, Uq < uq < Uq , and both Uq and Uq satisfy (|3.18p . Now we set 

wf{x,t) = y(^, 4{x)±e'Ce{e^'/'" - 1)). 

Then the same argument as in Lemma 13.81 shows that {w~ , wf ) is a pair of 
sub- and super-solutions for problem (-P^). Furthermore, since w~{x,0) = 
1^0 (x) < "^0(3;) < Uq{x) = ■wf{x,0), the comparison principle asserts that 

wj{x,t) <u''{x,t) <wf{x,t) for x G n, < t < ^"^e^|lne|. (3.21) 
3.3 Proof of Theorem D 

In order to prove Theorem l3.1l we first present basic estimates of the function 
Y after a time of order r ~ | lne|. 

Lemma 3.9. Let rj £ (0,r/o) be arbitrary; there exist positive constants 
and Cj such that, for all e G (0,eo), 

(i) for alHe{-2Co,2Co), 

a--r]< Y{i^r^\ In e\,^) <a+ + ??; (3.22) 

(ii) for all ^ G {—2Cq, 2Cq) such that \(, — a\ > Cje, we have that 

if e,>a + Cje then Y {ji~^\\ne\,i) > a+ - i], (3.23) 
if i<a-Cje then y(/i-^| lne|, ^ < a_ + ??. (3.24) 

Proof. We first prove ()3.23p . For i> a + Cje, as long as Y{t,£,) has not 
reached a+ — rj, we can use (|3.14p to deduce that 

Y{T,i) >a + Cie'^-(e-a) 

> a + CiCje^'^e 

>a+-r], 

provided that r satisfies 

1 , a+ — a — 1] 
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Choosing 

max(o — a_ , a+ — a) — rj 

= cl ' 

we see that /x~^|lne| > r^, which completes the proof of (|3.23|) . Using 
()3.15p . one easily proves ()3.24p . 

Next we prove (|3.22|) . First, in view of the profile of /, if we leave from a 
^ E [a_ — a++r/] then y(r,^) will remain in [a_ — r/, a+ +77]. Now suppose 
that a+ + ?/ < C < 2Co. We check below that Y{fi~^\ lne|,^) < a+ + ?/. First, 
in view of (jl.lj) , we can find p > such that 

if a+ <u< 2Co then /(n) < p{a+ - u) 

(3.25) 

if — 2Cq < u < a- then f{u) > —p{u — a_). 

We then use the ordinary differential equation to obtain, as long as a+ + 77 < 
Y < 2Co, the inequality Yr < — Y). It follows that 



< —p. 



Y-a^ 

Integrating this inequality from to r leads to 

Y{t,0 < a+ + - a+)e-P- 

< Q+ + {2Co - a+)e-P\ 

One easily checks that, for e £ (0, Eq)) with eq = eo(??) small enough, we 
have Y{ij,~^\ lne|,^) < q;_|_ + rj, which completes the proof of p.22p . □ 

We are now ready to prove Theorem 13.11 By setting t = fi~^e'^ \ lne| in 
(j3.2ip . we obtain 

y(/i-i| lne|,no (x) - (Cqe - Cee')) 

< u'{x,fi-^e'^\\ne\) < Y {n-^\lne\,u^ {x) + Cee - Cqe'^) . (3.26) 

Furthermore, by the definition of Co in (jl.Sp . we have, for eq small enough, 

-2Co < u^{x) ± {CqE - CqS^) < 2Co for x e n. 

Thus the assertion (13. 2p of Theorem 13.11 is a direct consequence of ()3.22p 
and 

Next we prove (13. 3p . We choose Mq large enough so that Mqe — CqE + 
Cqe'^ > CjE. Then, for any x G 17 such that Uq{x) > a + Mqe, we have 

Uq{x) - {Cqe - CaE^) >a + Mqe - Cqe + C^e'^ >a + C7E. 

Combining this, (|3.26p and (|3.23p . we see that 

u'^{x, h^^e'^I Inel) > q+ — rj, 

for any x £ Q with Uq (x) > a + Mqe. From the definition of Uq it is clear 
that Uq (x) > a + Mqe if and only if uo{x) > a + Mqe, provided that e is 
small enough. This proves ()3.3p . The inequality ()3.4p can be shown the 
same way. This completes the proof of Theorem 13.11 □ 
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3.4 Optimality of the generation time 

To conclude this section we show that the generation time := /i~^e^| lne| 
that appears in Theorem 13.11 is optimal. In other words, the interface will 
not be fully developed until t comes close to t^. 

Proposition 3.10. Denote by t^-^ the smallest time such that (|1.14|) holds 
for all t G [^min' -^]- ^^6^ there exists a constant b = b{C) such that 

t'min > /^"^e^(|lne| -6), 

for all e G (0, eo). 

Proof. For simplicity, we deal with the case where ()3.18p is valid. In that 
case, ()3.20p holds for all small e > 0. For each 6 > 0, we put 

t^h) := ii-^e^{\\ne\-b), 

and evaluate u'^{x,t^{b)) at a point x G Q.^ where dist{x,VQ) = Ce. Since 
uq = a on Tt and since |Vtio| < Co by (11.80 . we have 

no(x) < a + CoCe. (3.27) 

It follows from this and (j3.14p that 

wt{x,t^{b)) =Y{ji^^{\\n£\ - b),UQ{x) + eCae"^ - e'^C^ 

<a + Gael i°=l-*(no(x) + eCee"^ - e'^Cq - a) 

<a + C2e~^e-\CoCe + eCee"^) 

= a + C2e-\CQC + CQe-^). 

Now we choose b to be sufficiently large, so that 

a + C2e-\CQC + C^e-^) < a+ - t]. 

Then the above estimate and ()3.20p yield 

u%x,t%b)) < wf{x,t%b)) <a+-r]. 

This implies that (|1.14p does not hold at t = t'^{b), hence t^{b) < The 
lemma is proved. □ 

4 Generation of interface in the general case 

In this section we extend Theorem 13. II to the case where ^ 0. The proof is 
more technical than the case = 0, but the underlying ideas are the same. 
Hence we will basically follow the argument of Section [3l simply pointing 
out the main differences. 
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4.1 The perturbed ordinary differential equation 

We first consider a slightly perturbed nonlinearity: 

fs{u) = f{u) + 6, 

where 6 is any constant. For \6\ small enough, this function is still bistable. 
More precisely, fs has the following properties, whose proof is omitted: 

Lemma 4.1. Let 6q be small enough. Then for any 6 £ {—5o,6o), 

(i) fs has exactly three zeros, namely a^{5) < a{6) < a^{5), and there 
exists a positive constant C such that 

\a-.{5) - a„| + \a{6) - a\ + \a+{6) - a+\ < C\6\. (4.1) 

(ii) We have 

fs>0 in (— oo, a_(5)) U (a((5), a+((5)), 
fsKO in {a^{S),a{5)) D {a+{5),+oo). 

(iii) There exists a positive constant, denoted again by C, such that 

H6)-f,\<C\6\, (4.3) 

where 

H{6) := fM^)) = na{6)). 

Now for each 6 G {—6o,6o), we define Y{t,S^;6) as the solution of the 
following ordinary differential equation: 

f y,(r, 6) = fs{Y{T, 6)) for r > 



(4.2) 



where ^ varies in (— 2Co,2Co), with Cq being the constant defined in (jl.Sp . 

To prove Theorem 13. H we will construct a pair of sub- and super- 
solutions for (P^) by simply replacing the function Y(t,^) in (j3.19p by 
Y{t,£,;6), with an appropriate choice of 6. For this strategy to work, we 
have to check that the basic properties of Y(t, in Subsection 13.11 carry 
over to Y{t,^;6). 

First, it is clear that all the differential and integral identities in Subsec- 
tion [XT] that follow directly from (|3.5p are still valid for (|4.4p . In particular, 
Lemmas 13.21 and 13.31 remain to hold if we replace Y{t,(,) hy Y{t,^] 6), f by 
fs and A{t,C) by A{t,^;6), where 

Next let us show that the basic estimates which we have established in 
Subsection 13.11 are also valid for Y{t,S^;5). The following lemma, which is 
an analogue of Lemma 13.41 is fundamental. 
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Lemma 4.2. Letr] G (0,?7o) be arbitrary. Then there exist positive constants 
60 = 6o{r]), Ci = Ci{r]), C2 = 02(11) o,nd C3 = C^ij]) such that, for any 
6 G [-60,60], 

(i) if ^ ^ ~ ^) then, for every r > such that Y{t, ^; 5) remains 
in the interval {a{6), a+ — rf), we have 

C-^eMWr < Y^{T,i-6) < C2e^^^^\ (4.5) 
|A(r,e;<5)| <C3(e'^('^)^-l); (4.6) 

(ii) the same estimates as above hold if the interval {a{6), a+ — rj) is re- 
placed by (a_ + rj, a{5)) . 

Proof. In view of (14. ip . we can choose a small constant 60 = 6o{r]) > 
such that (a((5),a+ — rj) C {a{6),a+{6)), for every 6 £ [—60,60]. Therefore 
fs{q) does not change sign in the interval (a(5),a+ — rj). Thus, in order 
to prove the lemma, we just have to write again the proof of Lemma 13. 4^ 
simply replacing Y{t,(,) by Y(t,^;6). We do not repeat the entire proof 
here. Instead, let us explain why Ci, C2 and C3 can be chosen independent 
of 6. In view of the proof of Lemma 13.41 it is sufficient to estimate, for 
q G ia{6), a+ — rj], the modulus of the quantity 

. , > f'{<i) - f'Mi)) 

by a constant depending on r], but not on 6 £ [-60,60]. Since 

we see that the function {q,6) ^ hs{q) is continuous in the compact region 
{ l^^l < 60, a{6) < g < a+ — r/}. It follows that |/i5(q')| is bounded as {q,6) 
varies in this region. This completes the proof of Lemma 14.21 □ 

Corollary 4.3. Let 1] G (0, ?7o) be arbitrary. Then there exist positive 
constants 60 = 6o{rj), Ci = Ci(ry) and C2 = C2{rj) such that, for any 
6 G [-60,60], 

(i) if i £ (a((5), a+ — rj) then, for every r > such that Y[t, ^; 6) remains 
in the interval (a((^),a+ — rj), we have 

Cie'^(^)-(e - a{5)) < Y{t,C;6) - a{6) < C2e^^^''>^C - am (4.7) 

(ii) if i £ (a- + rj, a{6)) then, for every r > such that Y{t, ^; 6) remains 
in the interval (a_ -\-rj,a{6)), we have 

C2e'^(^)"(e - a{5)) < Y{t, 6) - a{6) < de^^^^^^^ - a{6)). (4.8) 
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Proof. We can simply follow the proof of Corollary 13.51 In order to prove 
that Ci and C2 are independent of 5, all we have to do is to find constants 
Bi = Bi{rf) > and B2 = B2{'n) > such that, for all 5 G [-^o, ^0] and all 
q G {a{S),a+ - rj), 

Bi{q - a{d)) < fs{q) < B^iq - a{5)). (4.9) 

This can be easily done, since {q,5) 1— > f&{q)/{q — cl{5)) is a positive contin- 
uous function on the compact region { \5\ < 60, a{6) < q < a+ — □ 

Now, it is no trouble to establish an analogue of Lemmas 13.61 and 13.71 
with constants independent of 6. We claim, without proof, that: 

Lemma 4.4. Let ij G (0,r7o) and M > be arbitrary. Then there exist 
positive constants = ^oiVj^^) '^'^^ ^4 = C4{rj,M) such that, for any 
6 G [So, 60], 

(i) if ^ ^ [oi+ — rj, a+ + M], then, for all r > 0, Y{t, ^; 5) remains in the 
interval [a+ — r/, a+ + M] and 

\A{T,i;5)\ < C^T for t>0; (4.10) 

(ii) if ^ € [a- — M, a- + r]], then, for all t > 0, Y{t, 6) remains in the 
interval [a- — M, a_ + rj] and ()4.10p holds. 

Lemma 4.5. Let rj G (0, ?yo) be arbitrary and let Cq be the constant defined 
in (jl.Sp . Then there exist positive constants 60 = So{rj), C5 = C^{rj) such 
that, for all 6 G [—60, 60], for all t > and all ^ G (— 2Co, 2Co), 

|A(r,e;'5)| <C5(e^(^)--l). 

4.2 Construction of sub- and super-solutions 

We now construct a pair of sub- and super-solutions by modifying the defi- 
nition ([XT^ . We set 

wf{x, t) = y(^-^,uo{x) ± e^r{±eg, ^); ie^) 

where the function r(5, r) is given by 

r(5,r)=C6(e^(^)^-l), 
and the constant Q is chosen such that, for all small e > 0, 

|/(x,t,u)|<a for (x, t, u) G n X [0, T] X M, 
which, in view of (jl.Sp . is clearly possible. 
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Lemma 4.6. There exist positive constants and Cq such that for all 
£ G (0,eo); i'^s^jf^t) ^-^ ^ P^**^ of sub- and super- solutions for problem (P^), 
in the domain $7 x [0, /i~"'^e^| lne|], satisfying w~{x,0) = w~^{x,0) = uo{x). 

Proof. First, the same cut-off argument as in Subsection 3.2 enabies us 
to assume (I3.18P for simplicity. Hence satisfy the Neumann boundary 
conditions. We define an operator C by 

Cu := ut- Au- e^^(/(ii) - /(x, t, u)), 

and prove below that Cwf > by slightly modifying the argument which 
we have used to prove Cqw^ > in Section [3l A straightforward calculation 
yields 

1 



Yr-f{Y)+eg%x,t,Y) C^ix{eG)e 



t Yfc 
7^-Auo-^\Vuo\ 



If £o is sufficiently small, we note that ±eG S {—6o,So) and that, in the 
range < t < /i~"^e^| lne|, 

|e2c7e(e^'(±^5)*/^' - 1)1 < e2Cg(e"M(±ee)/M _ i) < c^, 
which implies that 

uo{x) ± eV(±eg, 4) G {-2Co, 2Co). 

These observations allow us to use the results of the previous subsection with 
the choices r := t/e'^, ^ := uq{x) + e'^r{eQ,t/£^) and 5 := eQ. In particular, 
the ordinary differential equation p.4p yields Yr = f{Y) +eQ, which implies 
that 



G + g'{x,t,Y) 



+ Y, 



CQfi{£Q)e 



Yf 



\Vuo\' 



By the choice of Q the first term of the right-hand side member is positive. 
Using the estimate oi A = Y^^/Y^ in Lemma 14.51 we obtain, for a constant 
C5 that is independent of e. 



Cw+ >Y^ 



C6/i(eg)e''(^^)*/=' - lAnol - C^ief'^'^We^ - l)|Vno|^ 



|A«o| +C5|Vnop 



In view of (14.30 . this inequality implies that, for e E (0,eo), with eq small 
enough, and for large enough, 



CwJ > 



^^2^ ~ CsC'o^ - Co 



> 0. 



This completes the proof of the lemma. □ 
Hence, as in Section^ the comparison principle can be applied to deduce 
w~{x,t) <u''{x,t) <wt{x,t) for 2; G < t < ^"^e^l lne|. (4.11) 
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4.3 Proof of Theorem 13.11 for the general case 

As in Subsection 3.3, we first present a key estimate of tlie function Y after 
a time interval of order r ~ | lne|. Roughly speaking, a perturbation 6 of 
order e does not affect the result of Lemma [ 



Lemma 4.7. Letr] G (0,?7o) be arbitrary. Then there exist positive constants 
£q and Cj such that, for all e G (0, Eq); 

(i) for alHe{-2Co,2Co), 

a--7]< Y{fi-'^\\ne\,^;±eg) <a+ + 7], (4.12) 

(ii) for all ^ G (— 2Co,2Co) such that \(, — a\ > Cje, we have that 

if iya + C-je then lne|, ^; ieC?) > a+ - r/, (4.13) 

if i<a-Cj£ then y(^^i| lne|,^; zbe^) < a_ + ??. (4.14) 

Proof. In the sequel, by e we always mean e G (0, eo)) with eq = eo{r]) 
small enough. In view of ()4.ip . we have, for C7 large enough, a + Cje > 
a{±eg) + ^CjE. Hence for > a + Cje, as long as Y{T,(,;±eQ) has not 
reached a+ — t], we can use (|4.7|) to deduce, as in Section O that (j4.13p is 
valid provided that 

T > In J—— =:fi ' e Ine , 

where mo = max(a — a_, 0+ — a). To complete the proof of (j4.13p we must 
choose C7 so that lne| — /i~^(e)| lne| > 0. A simple computation shows 
that 

n Ine -u (e) Ine = ; — Ine ; —In ^ . 



The first term, thanks to ()4.3p . is of order e| In e| . Hence, for C7 large enough, 
the upper quantity can be made positive for all e. The proof of (j4.14p is 
similar and omitted. 

Next we prove (j4.12p . First, we can assume that the stable zeros of 
f±eg, a-{-^eQ) and a+(±e^), are in [a- — r/, a+ + r/]. Hence, in view of the 
profile of /±eg, if we leave from a ^ G [a_ — f?, a+ + ?]] then y(r,,^; ite^) 
will remain in [a_ — 77, a+ + rj\. Now suppose that a+ + r/ < C ^ 2Co. We 
check below that y(/i~^| lne|, ^; ±e^) < a+ + As in Section [3l as long as 
a+ + r/ < y < 2Co, ([T^ leads to the inequality Yr < p{a+ -Y) + eG. It 
follows that 

< -p + e-, 

y — q:_)_ t] 
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which imphes, by integration from to r, that 

y(r, e; ±eQ) < a+ + (2C7o - a+)e^-^^'^^\ 

One easily checks that, for e, we have Y{pL^'^\\n£\,£^]^eQ) < a+ +r], which 
completes the proof of ()4.12p . □ 

We are now ready to prove Theorem 13.11 in the general case. By setting 
t = lne| in ()4.1ip . we get 

Y(^fi~^\ lne|,uo(x) - e'^r{-eQ , fj.~^\ lne|); -eG^ 

< u%x,fi-'^£^\lne\) < Y(^n~^\ln£\,uo{x) +e^r{eg,fi~^\ln£\);+egy 

(4.15) 

The point will be that, in view of (j4.3|) . 

^.^/z-M±££)ine = 0. (4.16) 

£—0 n 

It follows that 

62r(±eg,^-i|lne|) =C76e(6('^-^(±^^»/^ -e) G (^Cge, ^Qe). 

Hence, as in Section[3l the result (j3.2p of Theorem l3.1l is a direct consequence 
of (liT2]) and (fiT5]) . 

Next we prove (j3.3p . We take x G J] such that no(x) > a + Mqe; then 

tio(x) -e2r(-eg,/i-i(e)|lne|) > a + Moe - fCee 

> a + Cye, 

if we choose Mq large enough. Using ()4.15p and ()4.13p we obtain (13. 3p which 
completes the proof of Theorem 13. 1[ □ 



5 Motion of interface 

In Sections [3] and HI we have proved that the solution develops a clear 
transition layer within a very short time. The aim of the present section is 
to show that, once such a clear transition layer is formed, it persists for the 
rest of time and that its law of motion is well approximated by the interface 
equation (P^). 

Let us formulate the above assertion more clearly. By taking the first 
two terms of the formal asymptotic expansion ()2.3p , we get a formal approx- 
imation of a solution up to order e : 

u'{x,t) « ff{x,t):=Uo(^^^)+eUi(^x,t,^^y (5.1) 
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Here Uq, Ui are as defined in (12. 5p and ()2.15p . The right-hand side has a 
clear transition layer which lies exactly on Tt- Our goal is to show that this 
function is a good approximation of a real solution; more precisely: 

If becomes close to -uF at some t = to, then it stays close to 
for the rest of time. Consequently, Tf evolves roughly like Fj. 

In order to prove this assertion, we will construct a pair of sub- and 
super-solutions u~ and for problem (-P^) by slightly modifying the above 
function uF . It then follows that, if the solution satisfies 

u~{x,tQ) < u''{x,to) < uf{x,to), 

for some tQ > 0, then 

u~{x,t) < u^{x,t) < uf{x,t), 

for tQ < t <T, which implies that the solution stays close to u^. 

The rest of this section is devoted to the construction of these sub- and 
super-solutions. We begin with some preparations. 

5.1 A modified signed distance function 

For our later analysis, it is convenient to introduce a "cut-off signed distance 
function" d, which is defined as follows. First, choose do > small enough so 
that the signed distance function d defined in ()2.ip is smooth in the following 
tubular neighborhood of F: 

{(x,t) G Q^, \d{x,t)\ <3do}, 

and that 

distiVt, on) > Mo for all t G [0, T]. (5.2) 
Next let C(s) be a smooth increasing function on M such that 

s if \s\ < do 

C(s) = < -2do if s < -2do 
2do if s > 2do. 

We then define the cut-off signed distance function d by 

d{x,t) = C{dix,t)). (5.3) 

Note that |V(i| = 1 in the region {{x,t) G Qt, \d{x,t)\ < do} and that, in 
view of (15. 2p . Vd = in a neighborhood of dQ. Note also that the equation 
of motion {P^), which is equivalent to (j2.13p . is now written as 

dt = Ad--f{x,t) on Ff, (5.4) 

where 'y{x,t) is the function defined in (j2.16p . 
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5.2 Construction of sub- and super-solutions 

As we stated earlier, we now construct sub- and super-solutions by modifying 
the function in ()5.ip . Concerning the second term Ui, which is defined 
in (|2.15|) . the terms AUi and Uu do not make sense as we only assume 
that g{-, -ju) S C^'^'^'^^ . In order to cope with this lack of smoothness, we 
replace C/i by a smooth function Uf, which is defined by 

UL, + f'{Uo{z))Uf = g'{x, t, Uo{z)) - Y{x, t)Uo'{z), 

(5.5) 

Uf{x,t, 0)=0, Ul{x,t,-) L°^i^), 

where 

7^(a;, t) = co(G^(x, t, a+) - G^ix, t, a_)), (5.6) 
with G^{x, t, s) = g^{x, t, r)dr. Thus C/f (a;, t, z) is a solution of (j2.9p with 

A = Al{x,t,z) ■.= g'{x,t,Uoiz))--f'ix,t)Uo'{z), (5.7) 

where the variables x, t, e are considered parameters. Using (jl.5p and the 
same arguments as in the end of Section [21 we obtain estimates analogous 
to (j2.18p and (j2.19p . with a constant M independent of e: 

\Ul{x,t,z)\ < M, \V^Uf{x,t,z)\ < M. (5.8) 

Moreover, g^ being in x and in t, AxUf and Uf^ are solutions of (j2.9p 
with A = AxAq and A = A^^, respectively. Thus, in view of ()1.3p . we obtain 

\AxU!{x,t,z)\ < C/e, \Uft{x,t,z)\ < C/e, (5.9) 

with some constant C independent of e. Similarly, (jl.5p and Lemma [2.3l vield 
estimates analogous to (I2.23P and ()2.24p for Uf, with C and M independent 
of e: 

M,{x,t,z)\ + \Ul,{x,t,z)\ < Ce-^l^l, (5.10) 

|V,C/f,(x,f,z)| <M. (5.11) 

In the rest of this section, C and M will stand for the constants that ap- 
pear in inequalities ()5.8p - ()5.1ip . Note also that (II. 7p implies the Neumann 
boundary conditions (|2.26p for Uf. 

We look for a pair of sub- and super-solutions uf for (P^) of the form 

nf{x,t) = Uo ( f ^^^^^ ) + eU! {x, t, f '^^'^ ) ± g(t) , (5.12) 

where 

g(t) = cT(/?e-^*/^' +e2Le^*). 
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Note that q = ae^ pt- It is clear from the definition of uf that 

lim.,±(x,t) = |"^ foran(x,t)eQ^ 
<^-*o 1 a- for all G Q^. 

The main result of this section is the following: 

Lemma 5.1. Choose /?, a > appropriately. Then for any K > 1, there ex- 
ist constants eq, L > such that, for any e £ (0, eo); the functions {u~,uf) 
are a pair of sub- and super-solutions for (P^) in the domain Q x [0, T]. 

5.3 Proof of lemma 15.11 

By virtue of p. 260 and the fact that Vd = near dQ, we have 

^ = on 9!^ X [0,r]. 

What we have to show is 

Cut ■■= «)t - A< - - eg'{x,t,ut)) > 0, 

and that Cu~ < 0. We will prove only the former inequality for uf, since 
the latter follows by the same argument. 

5.3.1 Computation of Cuf 

Straightforward computations yield 

{u+)t = Uo'ij + pt) + eUlt + Ulidt + ept) + qt 
Vu+ = Uo^ + eVUf + Ul.Vd 

Aut = Uo"^-^ + Uo'^ + eAU! + 2VC/f, • Vd + UfJ-^ + Uf^Ad, 

where the function Uq, as well as its derivatives, are evaluated at 2; = 
[d(x,t) + ep{t)^/e, whereas the function C/f, as well as its derivatives, are 

evaluated at (^x,t, (^d{x,t) + £p{t))/e^. Note that V and A stand for 
and A^, respectively. We also have 

/(n+) = /(C/o) + {eU! + q)nUo) + \{eUl + qf f" {6) 
g{x, t, uf) = g{x, t, Uo) + {ellf + q)gu{x, t, w), 
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where 6{x,t) and uj{x,t) are some functions satisfying Uq < 9 < uf , Uq < 
UJ < uf. Writing = g + — g and combining the above expressions with 
(12.51) and (15.51). we obtain 



Cut = Ei + --- + ^7, 

where: 

Ei = -^q [f'm + ^qf"{0)) + Uo'pt + Qt 

i^2 = (^ + ^)(l-|Vdp) 

E3= [^ + Ul){dt-Ad + j) 

= eUl.pt + p{gu{x,t,u;) - U!f"{e)) 

E5 = -7 Ul - ^ {Ulff"{e) + Ulgu{x, t,u;)-2 VUl ■ Vd 
Ee = eUlt - eAUl 

Ej = \{g' - g){x, t,ut)- - g){x, t,Uo) + ^(7" - 7)ix, t)Uo' . 

Before starting to estimate each of the above terms, let us present some 
useful inequalities. First, by assumption (jl.ip . there exist positive constants 
b, m such that 

f'{Uo{z)) < -m if Uoiz) G a_ + 6] U [a+ - b, a+]. (5.14) 

On the other hand, since the region {z G M | Uo{z) £ [a_ + 6, 0+ — 6] } is 
compact and since Uq' > on M, there exists a constant ai > such that 

Uo'{z) > ai if Uoiz) £ [a- +b,a+-b]. (5.15) 

We set 

/?=-, (5.16) 

and choose cr that satisfies 

< a < min (uo, cJi, 0-2), (5-17) 

where 

ai 1 4/3 



m + Fi' (3+V F2(/3 + l)' 

-^1 '■= ll/'llL°°(a_,a+)) E2 '■= ||/"||L°°(Q_-2,a++2)- 
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Combining (I5.14p and (IS.lSp . and considering that c < (Tq, we obtain 

Uoiz) - af'iUoiz)) >am for - oo < z < oo. (5.18) 

Now let i^T > 1 be arbitrary. In what follows we will show that Cuf > 
provided that the constants £q and L are appropriately chosen. We recall 
that a- < Uq < a+. We go on under the following assumption 

EqM < 1, egLe^^ < 1 . (5.19) 

Then, given any e £ (0,eo), we have e\Ui{x,t, z)\ < 1 and, since a < ai, 
< q{t) < 1, so that 



a--2<uf{x,t)<a+ + 2. (5.20) 



5.3.2 The term E- 



Direct computation gives 

= 4e-^*/^'(/-a/3) + Le^*(/ + eVL), 

where 

/ = C/^ - af'm - y r (^)(/?e-/^*/^' + e'le^'). 
In virtue of (|5.18p and (|5.20p . we have 

2 

I>cjm-^F2{l3 + e'^Le^^). 

Combining this, (I5.19P and the inequality cr < cj2, we obtain / > 2a(3. 
Consequently, we have 

5.3.3 The term E^ 

First, in the region where \d\ < do, we have |Vd| = 1, hence E2 = 0. Next 
we consider the region where \d\ > do. We deduce from Lemma [2. II and from 
that : 

lE^l < C(4 + -)e-^\''+^P\/' < ^e-A{rfo/.-H). 
We remark that 0<i^ — l<p< e^'^ + K. Consequently, if we assume 



e^^+K<-^, (5.21) 
2eo 



then — — Inl > — , so that 
e - 2e' 



\E2\ < ^e-^'^«/(2^) < C2 :-- 
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5.3.4 The term ^3 

By dEl) and (f2l6]) . we have 

{dt- Ad + -f){x,t) =0 on Ft = {x £ n, d{x,t) = 0}. 

Since 7 is of class C^~^^'^~ by virtue of (jl.Sp . we see that the interface Tt is 
of class C^^'''~2~. Therefore both Ad and are Lipschitz continuous near 
r^. It follows that there exists a constant > such that: 

\{dt - Ad + -f){x,t)\ < N\d{x,t)\ for all {x,t) G Qt- 

Applying Lemma 12.11 and the estimate (jS.lUp we deduce that 

lE^l < 2Arc7Me-A|d/e+p| 
e 

< 2iV(:7max5gK|^|e-^l«+Pl 

< 2NCmax{\p\,h. 

A 

Thus, recalling that \p\ < e^^ + we obtain 

\E^\ < Csie'^' + K) + C3' , 
where C3 := 2NC and C3' := 2iVC/A. 

5.3.5 The term 

In view of (II. 4p and (IS.lOp . both and are bounded by some constant 
C. Hence, substituting the expression for pt and q, we obtain 

|i?4| <C4(-/3e-^*/^'+eLe^*), 

where C4 := C + cr(C + AfFa). 

5.3.6 The term E5 

In view of (I2.16|) . the term I7I is bounded by co(a+ — a^)C on 17 x [0, T]. 
Using (|1.4p and (|5.1ip . we easily obtain |ii^5| < C5, where C5 depends only 
on C, M, F2. 

5.3.7 The term Eq 

We use dSj]) to deduce that \Eg\ < 2C =: Cq. 
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5.3.8 Finally the term Ej 

We recall that \g^ — g\ < Ce so that — 7I < co(a+ — a_)Ce. It then 
follows that 

\Ej\ <2C + Cco{a+ - a_) =: C7. 

5.3.9 Completion of the proof 

Collecting all these estimates gives 

Cut > - ^)e-^*/^' + {2a(3L - C3 - eCiL)e^' - Cg, (5.22) 

where Cg := C2 + KC3 + C3' + C5 + Ce + C7. Now we set 

T 4eo' 

which, for small enough, validates assumptions (j5.19p and (j5.2ip . For 
£0 small enough, the first term of the right-hand side of (I5.22p is positive, 
hence 

Cuf > [a(3L - C3]e^* - Cg > ^a(3L - Cg > 0. 

The proof of Lemma 15.11 is now complete, with the choice of the constants 
/?, a as in (f5J6]) . dSH]). □ 

6 Proof of the main results 
6.1 Proof of Theorem [TTSl 

Let 7] e (0,r/o) be arbitrary. Choose /3 and a that satisfy (15.160 . (15.170 and 

aP < |. (6.1) 



By Theorem 13.11 there exist positive constants £q and Mq such that (j3.2p . 
()3.3p and (13. 4p hold with the constant r/ replaced by a (5/2. Since Vuo • n ^ 
everywhere on Fq = {x G fi, uq{x) = a} and since Fq is a compact 
hypersurface, we can find a positive constant Mi such that 

if do{x)> MiE then uo(a;) > a + MqE 

(6.2) 

if (^0(2;) < — Mie then mo(x) < o — Mqe. 

Here (io(2;) := c?(x,0) denotes the signed distance function associated with 
the hypersurface Fq. Now we define functions H~^{x),H~{x) by 



H+{x) 
H^{x) 



a+ + a(5/2 if do(x) > -Mis 

a_ + a(3/2 if 4(2;) < -Mie, 

a+ - 0-/3/2 if dQ{x) > Mis 

a_ - 0-/3/2 if do(rE) < Mis. 
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Then from the above observation we see that 

H-{x) < u'{x,fi-^e'^\lne\) < for x £ Q. (6.3) 

Next we fix a sufficiently large constant K > such that 

Uoi-Mi +K)>a+-^ and Uo{Mi -K)<a^ + ^. (6.4) 



For this K, we choose Eq and L as in Lemma 15. II We claim that 

u-{x,0) < H-{x), i7+(x) < ti+(x,0) for X GO. (6.5) 

We only prove the former inequality, as the proof of the latter is virtually 
the same. Then it amounts to showing that 



[X . 



u;ix,0) = Uo{^^-K)+eUl{x,0,^^-K)-aif3 + e^L) < H 

(6.6) 

By (15. 8p we have l^/f | < M. Therefore, by choosing eq small enough so that 
EqM < a (3/6, we see that 

u~ {x, 0) < Uo{^^ -K)+eM- a{f3 + e^L) 
^ e ^6 

In the range where do{x) < Mie, the second inequality in ()6.4p and the fact 
that Uq is an increasing function imply 

Uq[ K)--al3 < - — = H {x). 

e 6 2 

On the other hand, in the range where do(x) > Mie, we have 

Uo{^ -K)- la(3 < a+ - \aP < H~{x). 
£ 6 b 



This proves (|6.6p . hence (j6.5p is established. 
Combining ()6.3p and ()6.5p . we obtain 



u~{x,{)) < u''{x,n~^e^\\ne\) < u+(x,0). 

Since u~ and uf are sub- and super-solutions of (P^) thanks to Lemma f5.lt 
the comparison principle yields 

u^{x,t) <u%x,t + t'') <uf{x,t) forO<^<^-^^ (6.7) 

where = /i~^e^| lne|. Note that, in view of (|5.13p . this is enough to prove 
Corollary 11.41 Now let C be a positive constant such that 

Uo{C - e^^ - K) > a+ - I and Uo{-C + e^^ + K) < a_ + |. (6.8) 
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One then easily checks, using ()6.7p and ()6.ip . that, for eq small enough, for 
< t < T - ^^ we have 



if d{x,t) > Ce then it^(x, t + t^) > 0+ — ?? 
if d{x,t)<-Ce then t + t^) < q_ + ry, 



(6.9) 



and 

u^{x,t + f) G [q_ - ?7,a+ + r/], 
which completes the proof of Theorem II. 3[ □ 

6.2 Proof of Theorem 11.51 

In the case where /i~^e^|lne| < t < T, the assertion of the theorem is a 
direct consequence of Theorem II. 31 Thus, all we have to consider is the case 
where < t < lne|. We first need the following lemma concerning 

Y, the solution of the perturbed ordinary differential equation (j4.4|) . 

Lemma 6.1. There exists a constant Cg > such that 

if C>a + C8e then y(r,^;±e^)>a for < r < ^~^| lne|, 

if ^<a-Cse then Y{T,^;±eg) < a for < r < ^^^1 lne|. 

(6.10) 

Proof. We only prove the first inequality. In view of estimates (|4.7p and 
()4.ip . we obtain, for ^ > a + Cge, 

y(r, ±eg) > a{±eg) + Cie^^^'^'^^a + Cge - a{±eg)) 
>a-Cge + Cii-Cge + Cse) 
>a + e{CiC8-Cg{Ci + l)) 
> a, 

if we choose Cs large enough. □ 

Now we turn to the proof of Theorem 11.51 We first claim that there 
exists a positive constant M2 such that for all t G [0, fi~^e'^ \ lne|], 

rf cAAA/2.(ro). (6.11) 

To see this, we choose Mq large enough, so that Mq > Cg + 2Ce holds in 
addition to ([32]), (f^ and ([331). We then choose M2 > Mi, where Mi is 
as defined in (j6.2p . In view of this last condition, we see that if eo is small 
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enough and if do{x) > M2e, then for < t < /i e \ lne\, 
uo{x) - e^r{-eg, ^) > a + Mqe - e^Ce [e^(-^5)l i^^l/^* - l] 

> a + e [Mo - Qe^''"^^^"^)^/^ + sCq] 

>a + e{Mo - 2Cq) ( ^ thanks to (ITOD ) 

> a + CsS. 

This inequahty and Lemma l6.ll imply ■w~{x,t) > a, where w~ is the sub- 
solution defined in (|3.19|) . Consequently, by (j3.2U|) . 

u^{x,t) > a if do{x) > A/^e. 

In the case where do{x) < —M2e, similar arguments lead to u^{x,t) < a. 
This completes the proof of ()6.1ip . Note that we have proved that, for all 

u'{x,t)>a if xen+\MM2e0^o), 

u^{x,t)<a if X G \7VAf2£(^o)• 
Since Ft depends on t smoothly, there is a constant C > such that, for all 
t G [0,//-^e2|lne|], 

roCAA^,.|i„,|(r,), (6.13) 

and 

nt \ Mc.iTt) c n+ \ MkUTo), 

As a consequence of ()6.1ip and (16.131) we get 

which completes the proof of Theorem II. 5[ □ 

Proof of Corollary 11.61 In view of Theorem 11.51 and the definition of the 
Hausdorff distance, to prove this corollary we only need to show that 

TtCAfc'em) for 0<t<r, (6.15) 

for some constant C" > 0. To that purpose let C' be a constant satisfying 
C > max(C,C), where C is as in Theorem 11.31 and C as in (j6.14p . Choose 
t G [0, T] and xq G Fj arbitrarily and, n being the Euclidian normal vector 
exterior to Ff at point xq, define a pair of points: 

a;"*" := xq + C'en and x^ := xq — C'en. 
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Since C > C and since the curvature of Tt is uniformly bounded as t varies 
over [0, T], we see that 

x+ G \ ^Ce(Xt) and X- £ \ Mce{rt), 

if e is sufficiently small. Therefore, if t G [;U~^e^| lne|,T], then, by Theorem 
11.31 we have 

u%x~,t) <a<u%x+,t). (6.16) 

On the other hand, if t £ [0,ij,-^e'^\ lne\], then from (Kl2h . (f6l^ and the 
fact that C > C, we again obtain (pT^ . Thus (f^T^ holds for all t G [0, T]. 
Now, by the mean value theorem, we see that for each t £ [0, T] there exists 
a point xi on the line segment x"*"] such that u^{xi,t) = a. This implies 
xi £ r|. Furthermore we have d{xo,xi) < C'e, since xi lies on the line 
segment This proves (j6.15p . □ 

7 Application to reaction-diffusion systems 

In this section we discuss the singular limit of the reaction-diffusion system 
{RD^) and prove Theorems 1 1 . 1 11 [1 . 1 3 1 and their corollaries. Our strategy is 
to regard the first equation of (RD^) as a perturbed Allen-Cahn equation 
and apply what we have already proved for this equation. 

7.1 Preliminaries: global existence 

Before studying the singular limit of (RD^), we first show that the solution 
of this system exists globally for t >0, provided that e is sufficiently small. 
Recall that the system (RD^) is written in the form 

ut = Au + ^ {f{u) + e /i(u, v) + 0{e^)) , 
_ vt = DAv + h{u, v), 

where h{u,v) satisfies the hypothesis (H). The standard parabolic theory 
guarantees the existence of local solutions for {RD'^). In order to prove that 
the solution exists globally for t > 0, it suffices to show that the solution 
remains uniformly bounded. This will be done by using the well-known 
method of invariant rectangles. 

Given arbitrary mq, vq £ C{Q), we choose a constant L > such that 

f{-L) > > /(L), -L < no(x) < L for x G a (7.1) 

Such a constant L exists since f{u) > for n < a_ and f{u) < for u > a+. 
By hypothesis (H), we can choose a constant Mi satisfying 

Ml > ||vo||L°°(f^), 
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along with the condition (ll.lSj) . namely 



h{u,-Mi) >0> h{u,Mi) for \u\ < L. (7.2) 

Now we consider the rectangle 

7^ := { {u,v) G I \u\ < L, \v\ < Mi }. 

It follows from (|7.ip that, for all sufficiently small e > 0, 

f{-L,v)>0>f{L,v) for|^;|<Mi. (7.3) 

The inequalities (|7.2|) and (|7.3|) imply that the rectangle 7^ is a positively 
invariant region for the system of ordinary differential equations 

ut = ^f{u,v), 
vt = h{u,v), 

since the vector field (e~^/''(ti, v), h{u, v)) points inwards everywhere on the 
boundary of TZ. The maximum principle then implies that TZ is also posi- 
tively invariant for the system (RD^). Consequently, since {uq{x),vq{x)) G 
TZ for X G fi, we have 

{u{x,t),v{x,t)) £TZ for X G n, t > 0, 

so long as the solution is defined. This uniform bound then implies that the 
solution exists globally for t > 0. 

In the case of equations for which only nonnegative solutions are to be 
considered (see Remark ll.9|) . we can argue just similarly, by replacing TZ by 
the rectangle TZ+ := {{u,v) \ 0<u<L,0<v< Mi}. Summarizing, we 
have proved the following proposition: 

Proposition 7.1. Let {uq,vq) G C{^}) x C(0). In the case where the condi- 
tions of Remark \1.9\ apply, assume further thatuo,vo > 0. Then there exists 
eo > such that for any e G (0, £q), the solution of (RD^) exists globally for 
t >0 and is uniformly bounded. 

Remark 7.2. For the details of the method of invariant rectangles, we refer 
the reader to the book |28t Chapter 14, Corollary 14.8]. See also [12] and 
|11| . It should be noted that [22] makes a much earlier study of invariant 
rectangles for a finite-difference scheme for reaction diffusion systems. 
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7.2 Proof of the main results 



Now we turn to the reaction-diffusion system (RD^) and explain our strategy 
for proving Theorems II. IH [1.131 and their corollaries. 

In what follows, we fix the initial data (uq, vq) and denote by (n^, f^) 
the solution of the system {RD^). The solution of the associated moving 
boundary problem (RD^) will be denoted by {T,v). 

Given a function v{x,t) on Q x [0, +oo), we set 

g^[v]{x,t,u) := -fi{u,v{x,t)) - e f^{u,v{x,t)), 

(7.4) 

g[v]{x,t,u) := -fi{u,v{x,t)), 

where /i, f2 are as in (jl.l7|) . The first equation of (RD^) is then written 
in the form 

ut = Au + ^{f{u)-eg'[v]{x,t,u)), (7.5) 

so that u^{x,t) is the solution of (P^) with the choice of the perturbation 
term g^(x,t,u) = g^[v'^](x,t,u). On the other hand, the equation of surface 
motion in the limit problem (RD^) is written in the form 

Vn = -{N -1)k + co g[v]{x,t,r)dr on Tj, (7.6) 

J a- 

SO that r is the solution of (P^) with the choice g{x,t,u) = g[v]{x,t,u). 

Thus Theorems 11.111 11.131 and their corollaries will follow from what we 
have shown for the single equation (-P^). In order for Theorems ll.3l and ll.5l 
for (-P"^) to be applicable to the present reaction-diffusion system {RD^)^ all 
we have to do is to verify the conditions ()1.3p to (II. 6p . More precisely, we 
have to show that, for all small e > 0, 

|A,/[t;^](x,t,n)| < Ce^i and \dt g'[v%x,t,u)\ < Ce^i, 

\dug'[v%x,t,u)\<C, 

\\f¥]{-r.u)\\ 1+^ _ < C, 
iiy L ' > - 

\g^[v''\{x,t,u) - g[v]{x,t,u) \ <Ce. 

Since g[v\, g^[v\ are defined by (17. 4p and since /i, f2 are smooth, it suffices 
to prove the following estimates for some C > and for all small e > 0: 

\A^v^{x,t)\<Ce-^ and \dtv^ {x,t)\ < Ce-^ , (7.7) 

Ib'^ll 1+^ _ < C, (7.8) 

\v''{x,t) -v{x,t)\ < Ce. (7.9) 

The estimates (|7.7p and (|7.8p are elementary, but (j7.9p requires far more 
elaborate analysis. In this subsection we prove ()7.7p . ()7.8p and give an 
outline of the proof of (j7.9p . A full proof of (|7.9p will be given later. 
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Proof of dZS]) and ([77]). Since v" satisfies 

vf = DAv^ + h{u',v'') inl7x(0,r], (7.10) 
along with the Neumann boundary conditions, it can be expressed as 

v%x,t) = h + h, (7.11) 

where 

h := / G{x,y,t)vo{y)dy, 
Jn 

h -= j / G{x,y,t - s)h{u^{y,s),v''{y,s))dyds, 
Jo Jn 

with G{x, y, t) being the fundamental solution for equation vt = DAv under 
the Neumann boundary conditions. Since h{u^,v^) is uniformly bounded, 
standard estimates of G{x,y,t) imply (|7.8p for any ■& £ (0, 1). 

In the mean while, the same rescaling argument as in Remark 11.71 yields 

lln^ll . ^ - < Ce~^. (7.12) 

Indeed, since VyU, Ur are bounded, where y = x/e, r = t/e^, we have 
S/xU = 0(l/e), ut = 0(l/e^). Consequently we have 

\u{x,t) — u{x' ,t')\ \u{x , t) — u{x' , t)\ \u{x' ,t) — u{x' ,t')\ 



\x-x'\^ + \t-t'f/^ ~ Ix-x'l"^ \t-t'f/^ 

, , , , , Ml ,q\u(x,t) — U(x' ,t)\^ 

< u{x,t) - u{x',t)^-^^-^^ jrw^^ 

\x — x'\^ 

+ln(x^t)-n(x^tol^--^/^ '"^"'^|:;;|;V^^''^' 

< (2||u||ioo)i-''||V,«||i^ + (2||n||Loo)i-'^/2||^,,||^/i 

Combining (17121) and (EHl), we see that ||/i(^t^ t;^)|| ^ _ < Ce"'^, 

hence, by the Schauder estimate, 

I I /q I I 4 1 ^ "J ^ CE ^ . 

II ^ 11^2+^,1+7 (f^x[0,T]) - 

Here the constant G may depend on the choice of G (0, 1). On the other 
hand, Ii is bounded in C^'^(r2 x [0, T]) since vq S G'^{Q). Combining these, 
we obtain \ AxV'^{x,t)\ = 0{e^^), hence 0(e~-^). Substituting this into (j7.10p 
yields the second inequality in (17. 7p . □ 
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Outline of the proof of ()7.9p . We decouple the system {RD'^) as follows. 
As mentioned earlier, u^[x^t) is the solution of {P^) with the choice of the 
perturbation term g^{x,t,u) = g'^[v^]{x,t,u), that is, 



ut = Au + ^(/(w) - eg''{x,t,u)) 

p = o 

Oh' 



^u{x,0) = uo{x). 

Once the solution is determined, is the solution of the problem 



(^) 



vt = DAv + h{u, v) 

^ V{x,{)) = Vq{x), 



with the choice u = . This means that v'^{x,t) is a fixed point of the 
following map := ^2° ^i'- 



$f via (it) _ $2 via 
<P : V > U 



> V 



where $f maps a function v{x,t) to the solution u^(x,t) of (^) for the choice 
g'^{x,t,u) = g^[v]{x,t,u), and $2 maps a function u{x,t) to the solution 
v{x,t) of for the choice u = u. 

On the other hand, as for the limit problem (RD^), the solution v{x,t) 
can be regarded as a fixed point of the map := ^2 ° '■ 



Q '^'j' via ^ <I>2 via 

q? : V > u > V, 



where maps a function v{x,t) to the step function 
u(x, t) 

where Tt [v] is the solution of the equation of surface motion 

h 

fi{r,v{x,t))dr on Tt 



a+ in n^{Tt[v]) 
a- in n~{Tt[v]), 



li=0 



-(A^-l)«-co 

= ro, 



and 17 (r) denotes the region enclosed by the hypersurface T and i^'*'(r) 
the region between 9fi and T. 

In what follows we set 

Qt:=nx {0,t) for < t < T. 
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Given 60 > 0, we define tmax = tmax{So) > by 

tmax = max{t G [0,T], Wv" - w||L°°(Qt) < So}- (7.13) 
The key estimates for proving (j7.9p are the following: 

Claim 7.3. There exist constants (5o > and C > such that, for any 
t G {0,tmax], we have 

W^^ivn - '&°(i))||Loo(QO < Cl^-j2=\\v' - i\\Loo(^Q^)ds. (7.14) 

Claim 7.4. There exists a constant ^ > such that, for any v satisfying 
the estimates ()7.7p . (I7.8p and the Neumann boundary conditions, and for 
any t G (0, T], we have 

||$^(i;)-$°(^)|Uoo(Q,)<Ae. (7.15) 

The proof of these claims will be given later. For the moment, let us 
simply mention that Claim [77^ can be shown by the following two-step ar- 
gument: first, our results on the single equation (P^) yields 

<^Uv)-<^'i{v) = 0{e), 

in the sense that the transition layer of $i(f) and that of within 
an 0(e) distance; this observation and an estimate of the heat kernel yield 
()7.15p . To prove Claim 17.31 we also use a similar estimate of the heat kernel, 
see Lemma 17.61 and Subsection 17.71 for details. 

Combining these estimates, we obtain, for any t £ {0,tjnax], 

/■* 1 

< Ae + C / 11^;'' - v\\Loo,Q\ds. 

Jo Vt-s 

As we will see later in Lemma 17.71 this implies 

Ik" - v\\loo(q,) < Aek{t) (0 < t < Wx), (7.16) 
where k is the function determined by the integral equality: 

k{t) = l + C [' 4^ds. (7.17) 
Jo Vt-s 

Since k{t) is bounded on any finite interval [0, T], we obtain — "yllLoo^g^) = 
0{e), for < t < t 

max- This implies, first of all, that t^ax — T if e is small 
enough, hence it proves (j7.9p . for < t < T. □ 



The rest of this section gives a detailed account of the proof of ()7.9p . We 
begin with some notations to clarify the statements of the above claims. 
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7.3 Some notations 



Given any function g{x,t,u) satisfying the conditions ()1.3|) and (II. 5p . we 
can define a classical solution of the interface equation (P^) on some time 
interval < t < T"*"^'(^). We denote this solution by Tt[g] in order to clarify 
its dependence on g. More precisely, Tt[g] is a solution of the problem 



Vn 

Til 



i=0 



-{N -1)k + co 

= ro. 



g{x^ t, r)dr on Tt 



Also, we denote by u'^[g\{x,t) the solution of the problem 



{PI) 



( 1 

Ut = H — 7^{f{u) — £g{x, t, u)) in O x (0, +oo) 

— = on X (0, +oo) 



0) = uo{x) 



in 17. 



Once the interface Tt[g] is given, we denote by 17^ [^], 0,f[g] the region en- 
closed by Tt[g] and the one enclosed between dO, and Tt[g], respectively. As 
in (jl.lSp . we define the step function u[g]{x,t) by 



u[g]ix,t) 



a+ in 17^ 



for t G [0,T""'^g)). 



(7.18) 



Next, given any function u{x,t) on 17 x [0, T], we denote by 
the solution of the problem 



(7.19) 



Vt = DAV + h{u{x,t),V) inl7x(0,r] 
dV 

— = on ai7 X (0, T] 

dv 

V {x,{)) = vq{x) in 17. 

In view of (j7.4p and the above notations, the solution {u'',v^) of {RD 
is expressed as 

On the other hand the solution (r,?;) of (RD^) is expressed as 

rt = Tt[g[v]], v = V[u], 
the step function u in (RD^) being given by 

u = u[g[v]]. 

Finally, the maps <I>f , and <I>2 are now written as 

: V ^ u^g^v]], : V ^ u[g[v]], ^2 ■ u ^ V[u]. 



49 



7.4 Interface motion under perturbation 

In this subsection we show that the interface Tt[g] depends continuously on 
the pressure term induced by g. To this end, we first fix constants C^, > 0, 
i!) G (0, 1), and denote by y the set of functions g{x, t, on x [0, T*] x R 
satisfying 

sup 1+^ <C,. (7.20) 

Proposition 7.5. Let g € y. Let T £ (0, r"*"^(5)). Then there exist 
positive constants 6, K, M such that, for any g £ y satisfying 

\\g - 5llL°°(nx(o,T)xE) < (7.21) 

it holds that 7"™"^(^) > T, where we recall that x"^°'^(g) is the maximum 
time of existence of a classical solution of Problem (P^)- Furthermore, for 
each t G [0,r], 

dn{rt[g],rt[g] ) < K{e^' - 1) - 5llL-(Qx(o,t)xM) • (7.22) 

Proof. First, the assertion that T^'^^{g) > T, for g sufficiently close to g, 
follows from the standard local existence theory for quasi-linear parabolic 
equations. In fact, by using appropriate parametrization, one can express 
Tt[g] and r([^], as graphes over M, where M is a A^— 1 dimensional manifold 
without boundary, and transfer the motion equations (P^) and (P^), into 
quasi-linear parabolic equations on the manifold M, at least locally in time. 
For more details we refer to [9]. Since g and g satisfy ()7.20p . and since the 
embedding 

is compact if < < ??, the assumption ||^ — ^Hloo < 6 implies 

\\gi-,-,u) - g{-,-,u)\\ < c{6), 



where c{5) is a constant satisfying c{6) — > 0, as ^ 0. Consequently, the 

t) and (Pi* 



coefficients appearing in {P^) and {P^) satisfy 



g{-,-,r)dr- / g{-,-,r)dr \\ , ^+^, < {a+ - a-)c{5). 

J a- ^ ^ 

Hence, the two solutions Tt\g] and Tt\g\ stay close to each other, at least lo- 
cally in time, and, by repeating this argument, one can prove that 7"™-<^^(^) > 
T, for 5 sufficiently small. 

Next we prove the estimate ()7.22p . This will be done by using the maxi- 
mum principle. Let us introduce some notation. For each g G 3^, we denote 
by d{x,t;g) the signed distance function associated with the interface Tt[g]. 
By Tf < Tt we mean that F^ lies inside of F^. Clearly we have 

Tt [g] ^Ttig] ^ d{x,t;g)>d{x,t;g) for xeU. (7.23) 
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Now we choose to £ [0, T] arbitrarily and put 

Vo ■= \\g - 9\\L^(nx{0,to)xR)- 

Then 

g{x, t, u) - T]o < g{x,t,u) < g{x,t,u) +7J0 for 0<t<to. 

The comparison principle then yields 

rt[9-Vo]^rt[g]^Tt[g + 7]o] for < t < to- 

Thus, in order to prove (17.220 . it suffices to show that there exist constants 
K,M > such that, for all small tjq > 0, 

r dn{Tt[g-rio],Tt[g]) < Kr,o{e^''-l) 

{ (7.24) 

[ dniTtig + r/o], TM) < Krioie^"' - 1), 

for < t < to- We will only show the latter inequality for rt[^ + tjq] since 
the former can be shown in the same manner. 

Recall that d{x,t;g) satisfies the equation (|5.4p . namely 

dt = Ad-Co / g{x,t,r)dr on r([^]. (7.25) 

J a— 

Choose a constant do > such that d{x, t; g) is smooth — say, in x and 
in t — in the neighborhood A/^q (Fj [5] ) , < t < T. By the smoothness 
of d{x,t;g) and equality ()7.25p . there exists a constant N > such that 

\dt-Ad + co [ ^ g{x,t,r)dr\ < N\d\ in Md,/2{rt[g]). 

Now we put 

d"^"'(x,t) ■.= d{x,t;g)-Krjo{e^'''-l), 
ft :={x £n\d''^'"{x,t) =0}, 

where the constant K is to be determined later. If 

^ * e do 1 

m < Vo ■= — 4 — K , 

then Tt lies within the neighborhood J\fda/2(^t[g])- Observe that 
(d"^"')t - Ad"^"' = dt- 2NKrioe^^^ - Ad 

< -co / ^ g{x,t,r)dr + N\d\ -2NKr]oe'^^\ 
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Since d = Kr]o (e^^* — 1) on Tt, we obtain 

{d^'^)t - Ad^'^ < -co n g{x, t, r)dr - NKrjo on f*. 

J a- 

Now we set 

K = {a+- a_)coiV"^ 
Then it follows from the above inequality that 

(d"^"-)* < Ad"^'" - CO n g{x,t,r)dr - (a+ - Q-)cor?o on f*. 

J CL- 

This inequality and the fact that d"'^^{x, 0) = d{x, ; g) imply that Tt satis- 
fies 

Vn> -{N - 1)k + Co {g{x,t,r) +rio) dr on Tt, 

J a- 



\t=0 



To. 



On the other hand, Tt[g + 7?o] satisfies 

Vn = -{N - 1)k + Co {g{x,t,r) +r]o) dr onTtig + m], 

J a- 

_^t[g + vo]\t=o = ^0- 
By the comparison principle, we obtain 

Tt[g]^Tt[g + vo] ^ f f for < t < to- 

Consequently, 

dn{rt[g + vo], Ttig]) < dn{tt, Tt[g]) < K^o (e'"^* - 1), 
for < t < to- The lemma is proved. □ 

7.5 Proof of Claim YTM 

For a function v{x,t) satisfying the estimates (j7.7p . (j7.8p and the Neu- 
mann boundary conditions, we compare below ^^[v) = <I>2 ° ^li'^) and 

Action of <I>f and <I>5' Let us compare $f(t;) = with the step 

function <^\{v) = u[g[v\\. By the definitions in ()7.4p we have g'^[v\ = g[v] + 
0(e), and all the conditions in (jl.3p - (|1.7p are satisfied. It follows that our 
results for the single equation apply and, in particular, 

u-{x,t) < u''[g''[v]]{x,t + t^) < uf{x,t) foTO<t<T-f, 
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where are as in (15.120 . d being the signed distance function associated 
with the interface rt[g'[i;]]. Since the term e~^*/^ in q{t) — that appears in 
()5.12p — quickly becomes smaU, 

n^[/H] - u[g[v]]\{x,t) < a+ - u,(^ d{x,t) - ep{t) ^ ^ ^^^^ ^^^^^^^^ 
u'mv]]-u[g[v]]\i^,t) < Uo[^^^^^^^^)-a^+0{e) in n^[g[v]] 

for /iie^l lne| <t<T, provided that we choose the constant fii large enough. 
Consequently, by Lemma [2711 there exist constants B, C > such that 

\mv) -^^^{v)\{x,t) < Sexp(-Ai^^^) +Ce, (7.26) 

for {x,t) e n X [fiie^\lne\, T]. 

Action of $2. Next we compare ^'^{v) = V[^l{v)] and ^'^{v) = V[^i{v)]. 
Set w := ^''(v) — ^^{v). By subtracting the equations for y[<I>|(t;)] and 
we obtain 



wt 



DAw+(^h{mv),^'{v))-h{^'i{v),^°{v))y 



Since |/i($f(v),$^(?;))-/i($?(?;),^'°(t;))| < C\w\ + C\^l{v) -^°^{v)\ for some 
constant C > 0, the function w := e~^^w satisfies 

Wt < DAw + Ce^^^\^l{v){x,t) - ^'({v){x,t)\ +C{\w\ -w), 

hence 

Wt < DAw + C|$f (t;)(x, t) - $?(^;)(x, t)\ + C{\w\ - w). (7.27) 

Now let W{x,t) be the solution of the equation 

Wt = DAW + C\mv){x,t) - ^'({v){x,t)\ +C{\W\ -W), 

with initial data W{x,0) = 0. Then since (17.271) implies that is a sub- 
solution of the above equation, and since w{x,0) = 0, we have 

w{x,t) <W{x,t) forxeH, t>0. (7.28) 

Moreover, since W > 0, the above equation for W can be reduced to 

Wt = DAW + C\^l{v){x,t) - <^^^{v){x,t)\. 

In view of this, we see that 

W{x,t) = C f I G{x,y,t - s)m{v){y,s) - ^\{v){y,s)\dyds, 
Jo Jn 
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G{x,y,t) being the fundamental solution that appears in (17.111) . This and 
(|738D yield 



r-t 



\w{x,t)\<Ce'^' / G{x,y,t-sW,{v){y,s)-^^,{v){y,s)\dyds. (7.29) 
Jo Jn 

Combining this and ()7.26p . we obtain 

\w{x,t)\ < BCe^' f [ G{x, y, t-s) exp (-A ^^^ )dyds+0{e). (7.30) 
Jo Jn V e / 

In order to estimate the above integral, we need the following lemma: 

Lemma 7.6. Let T be a smooth closed hypersurface in and denote by d{x) 
the signed distance function associated with T. Then there exist constants 
C, ro > such that for any function r]{r) > on M, it holds that 

/ G(x, y, t) 7jid{y)) dy<^ r,{r) dr forO<t< T. (7.31) 

J\d\<ro yt J^ro 

The proof of this lemma will be given in the next subsection. As is 
easily seen from its proof, the above estimate remains to hold if F depends 
on t smoothly; in other words, the constant C can be chosen uniformly as 
r varies. Applying the above estimate to rt[g'[t>]], < t < T, we obtain 



/ G{x, y,t — s) exp ( — A 
Jn ^ 



\diy,s)\ 



e 



dy 

\d{y,s)\ 



/ + / G{x, y, t - s)exp( - X 

J\d\<ro J\d\>ro ^ 

e 



dy 



0( ' 



It follows from this and (|7.30p that 

\w{x,t)\=o(e [ / ds) +0{e) =0{e), (7.32) 
^ Jo Vt- s ' 

which completes the proof of Claim 17.41 □ 

7.6 Proof of Lemma 17.61 

We first show that 

/ G(x, y, t)dSy<-^ for X G 0, < t < T. (7.33) 
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It suffices to prove this estimate on a small interval [0,to]) since the esti- 
mate for the remaining interval [to, T] will follow by simply choosing a large 
constant C (since G is bounded for t large). Hereafter, we choose to suffi- 
ciently small. Then, for < t < to, G{x,y,t) is well approximated by the 
fundamental solution on the entire space M^: 

^=(4vrZ?V/2 ^"P(-^4^)- 

In particular, there exists a constant C > such that 

< G{x, y,t) <C Go{x, y, t) ior x,y eH, < t < to, 

(see, for example, [Ml Section IV.2]). Thus it suffices to prove (I7.33P for Go 
instead of G. 

Given x £ il., let xq be the point on F that is closest to x, and let n(xo) 
be the outward normal to T at xq. Then x — xq = d{x)n{xQ). Define 

Y:={yeR^ , yn(xo) = 0}, Yq = span{n{xo)) , 

where • denotes the Euclidean inner product in and span{w) the line 
spanned by the vector w. This gives an orthogonal decomposition = 
y © Yq) and xq + Y is the tangent hyperplane of F at xo. Since F is smooth, 
it is expressed locally as the graph of a map defined on a subset of Y. More 
precisely, there exist a smooth map h : Y ^ Yq and a constant 6 > such 
that h{0) = 0, V/i(0) = 0, and that 

S:={xo + y + h{y) , y G f, |y| < 5 } C F, 

(7.34) 

dist{xQ, T \ S) > 6. 
Now we decompose the integral ()7.33p for Go as 

/ ^' = JiJty^ ( /, + i\s ( " ^i^) ^^^) • 

Since |x — y| > \d{x)\ for every y G F and since 

\x - y\ > \ \x - xo\ - |y - xol I = I \d{x)\ -\y- xo\\, 

we have 

I I ^ \dix)\ + I \d{x)\ - |y - xol I ^ \y- xo| 



This and UiTMh yield 

\x -y\> - for y G F \ S. 
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Consequently, 

exp (- !ij^)dS, in, (7.35) 



Jr\s 



where |r| denotes the total area of T. 

On the other hand, for each y £ S, we can express y — xq as 

y-xo = y + h{y) {y e Y, h{y) e Yq), 

and Y can be identified with R^-^. Thus 

,~,M2 



/ \x - xq - y - h{y)\ ^ , \v7hf~\\2 

exp jT- VI + NHy)\ 



\y\<s ' 4Di 
Since V/i(0) = 0, there exists a constant Ci > such that 

\Vh{y)\ < Ci\y\ for \y\ < 6. (7.36) 
Note also that the orthogonality {x — xq — h{y)) _Ly implies 

\x-XQ-y - h{y)\^ = \x-xq- h{y)\'^ + \y\^ > \y\'^ . 
Combining these, we obtain 

/^exp(-^).., </ exp(-|^)Vl + CW# 



J\z\<^^t ^5 ^ 



dz, 



'\z\<Vi 

where z := y/\/t. Observe that, as t — > 0, 

/ e-\'\"l^''Jl + tCl\z\^dz ^ [ e-l^l'/4^dz = (4I)^)(^-^)/2. 

Consequently, 

^ /exor-^i^^ciS„<^ + ori.V 



Combining the estimate above and ()7.35p . we obtain 

Since F is a smooth compact hypersurface, its curvature is bounded. There- 
fore, the constants 5 and Ci that appear in (jT.SSp . (j7.36p can be chosen 
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independent of the choice of xq G F. Hence the above 0(1/ \/t) estimate 
is uniform with respect to the choice of x G 0. This proves the estimate 

Now, choose a sufficiently smah constant tq > such that the signed 
distance function d{x) is smooth in the region {d{x) < 2ro}. For each 
r G [— J^Oj'^o]) we define a hypersurface T{r) by 

r(r) := { x G fi, d{x) = r }. 

Then the curvatures of r(r) are uniformly bounded as r varies, which implies 
that there exists some constant C > such that 

/ Gix,y,t)dSy < for 0<t<T, r G [-ro,ro]. 
Jr{r) ■ Vt 

The estimate (I7.3ip now follows by integrating in r. □ 
7.7 Proof of Claim [TH] 

We compare below = $2 o ^i{v'^) and ^^{v) = ^2 ° ^i{v). 

Action of Let us compare the two step functions ^'({v'^) = u[g[v'']] 
and ^i{v) = u[g[v]]. We want to apply Proposition 17. 5t with g[v] and g[v^], 
playing the role of g and g, respectively. (Hence, the role of T"^"'^(g) is played 
by T"^""^ {g[v]) , which corresponds to T™'^^ in Lemma [l.lU|) . First, we choose 
> large enough so that both t, u) and t, u) satisfy (I7.20p . 

For T G (0,T'""^), we choose 6, K and M as in Proposition 17. 5i Next, we 
define Ki = max(„^^)g7^ \dyfi{u,v)\^ with TZ being the rectangle defined in 
Subsection 17.11 and 5q = 5/Ki. We observe that, using the definition of 
tmax in (|7.13p, 

\\9¥] - 9[v]\\L^(nx(a,tma.)^V) < Ki\\v^ - w||L-(Qt^„j < Ki5q = 6. 
By (I7.22p . it follows that, for any t G [0,tmax]5 

dniTt[g[v% Tt[g[v]])< K{e^'' " 1) hV] - g[v]\\L- ■ 
Combining these, we obtain 

dH{Tt[g[v']l Tt[g[v]] ) < KK^{e^'' - 1) ||^;^ - ^||l^(q,) . (7.37) 

Action of <1>2' Next we compare the two functions <I>'^(u^) = y[$5(^^)] = 
y[u[5[i;^]]] and ^^{v) = V[<^\{v)\ = V[u[g[v]\\. Since 

\u[g[v'']] - u[g[v]]\ < a+ - a., 



57 



and since the two step functions differ only in the region enclosed between 
the two surfaces r([(7[i;^]] and rt[(jr[-u]], the estimates ()7.29p and (|7.3ip imply 
that there exists a constant Bi > such that 

\\V[u[g[v ]]] - V[u[g[v\]\\\L^^Q^) < Bi ds. 

Combining this and (j7.37p . we obtain, for any t G [0, tmax]! 

\\^0{vn-^''{v)\\L^^Q,)<C f^=\\v^-vU^^Q^)ds, (7.38) 

Jo vt — s 

with C = BiKKiie^''^ - 1). The proof of ClaimOis complete. □ 
7.8 Estimate of k{t) 

In this subsection we justify the estimate (|7.16p . Let k{t) be the function 
satisfying ()7.17p . namely, 

k{t) = 1 + C [ ds for t > 0. 

Jo Vt- s 

We will show below that k is given by 

^(t) = e^"*(l + C / —^ds). (7.39) 



The following lemma justifies (|7.16p : 

Lemma 7.7. Let k{t) be a continuous function satisfying 

G<k{t)<A + cf-^^ds {0<t<T), 
Jo Vt- s 

for some constant A > and T > 0. Then 

< k{t) < Ak{t) for 0<t<T. (7.40) 

Proof. Define 

k^{t) ■= {l + e)Ak{t). 

Then this function satisfies 

k,{t) = {l + e)A + C f -^0= ds (0 < t < T). (7.41) 
Jo Vt- s 

In particular, we have kg^O) = (1 + e)A > A> k{0). Let us show that 

k{t) < ke{t) for < t < T. (7.42) 
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Suppose that (I7.42|) does not hold. Then there exists to G (0, T] such that 

k{t)<ke{t) forO<t<to, k{to) = ke{to)- (7.43) 
Combining the first part of fTiS]) and (17^ . we get 

hito) = il + e)A + C r >A + C r > fc(to), 

Jo V *o - s Jo V *o - s 

but this contradicts the second part of ()7.43p . establishing ()7.42p . Letting 
e ^ 0, we obtain (17301). □ 

Corollary 7.8. Let k{t) he a continuous function satisfying 

< k{t) =A + C f -^0= ds {0<t<T), 
Jo Vt- s 

for some constant ^ > and T > 0. Then k[t) = Ak{t). In particular, the 
function k[t) is uniquely determined by the integral identity (|7.17p . 

Proof. Define k{t) := A-'^k{t). Then k{t) satisfies 

k{t) = l + C [' ^0=ds {0<t<T). 
Jo Vt- s 

By Lemma 17.71 we have k{t) < k{t). Exchanging the role of k and k, we 
obtain the opposite inequality, hence k{t) = k{t). Thus k{t) = Ak{t). □ 

Now let us prove (I7.39p . Integration by parts gives 

k{t) = 1 + 2CVt + 2C I Vt^ k'{s)ds. 
Jo 

Hence 

F(t) = -^ + C f'lkLds. 
Vt Jo Vt-s 

Thus the function m{t) := k'{t) — CjVt satisfies 

Jo Vt- s Jo Vs{t - s) Jo Vt- s 

Since the first integral on the right-hand side is equal to vr, we obtain 

m{t) = C\ + C t ^^^^ds. 

Jo Vt-s 

It follows from Corollary 17.81 that m{t) = C^7rfc(t), hence 

k'{t) = CVfc(t) + 

Vi 

This and the equahty Ai(0) = 1 yield (1739]) . □ 
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